Microphysics of 
Cosmic Strings in 
Supersymmetric and 
Grand Unified Theories 



Stephen Christopher Davis 

Trinity College, Cambridge 



A dissertation submitted for the degree of Doctor of Philosophy 
University of Cambridge 
September 1998 



Acknowledgements 



I wish to thank my supervisor, Dr. Anne C. Davis, for all her advice and support. I am 
also grateful to Dr. Mark Trodden and Dr. Warren B. Perkins for useful discussions. 



Declaration 

Apart from the introductory sections, including chapter 1, and where indicated in the 
text, the work presented is my own. The material in chapters 2-5 is published as follows: 
Parts of chapters 2 and 3 appeared in Physical Review D 55 (1997) 1879-1895 (with Dr. 
A. C. Davis), most of chapter 3 appeared in Physics Letters B 408 (1997) 81-90 (with Dr. 
A. C. Davis and Dr. W. B. Perkins), chapter 4 was published in Physics Letters B 405 
(1997) 257-264 (with Dr. A. C. Davis and Dr. M. Trodden), and chapter 5 was published 
in Physical Review D 57 (1998) 5184-5188 (with Dr. A. C. Davis and Dr. M. Trodden). 
The work contained in chapter 6 is in a preliminary form. Whilst the papers have been 
published in collaboration, I have only included my own contributions in this dissertation. 

To the best of my knowledge, the work presented as my own is original. 



i 



Contents 



1 Introduction 1 

1.1 Particle Cosmology and Cosmic Strings 1 

1.2 Grand Unification and Symmetry Breaking 3 

1.3 Topological Defects 6 

1.4 The Abelian Cosmic String Model 7 

1.5 Determination of the String Field Profiles 9 

1.6 String Superconductivity 11 

1.7 Evolution of Cosmic String Networks 13 

1.8 Overview 15 

2 Microphysics of SO(10) Cosmic Strings 18 

2.1 Introduction 18 

2.2 An SO (10) Grand Unified Theory 19 

2.3 GUT Strings 21 

2.4 Cosmic Strings and Multiple Higgs Fields 25 

2.5 The Electroweak Symmetry Breaking 26 

2.6 Electroweak Cosmic String Solutions 30 

2.6.1 The Abelian U(l) String 30 

2.6.2 X s and X' SU(2) Strings 32 

2.6.3 W R SU(2) Strings 33 

2.6.4 Y' SU(2) Strings 35 

ii 



2.6.5 Summary 39 

2.7 Other Symmetry Restorations 40 

2.7.1 The Intermediate Symmetry Restoration 40 

2.7.2 The Minimal Energy Choice of T s 41 

2.7.3 Non-gauge field symmetry restoration 41 

2.8 Other Related Grand Unified Theories 42 

2.9 Conclusions 44 

3 A Fermion Zero Mode Index Theorem 47 

3.1 Introduction 47 

3.2 Zero Mode Index 48 

3.3 Fermion Zero Modes on SO(10) Cosmic Strings 55 

3.3.1 Zero Modes for the Abelian String 56 

3.3.2 Zero Modes of the SU(2) Strings 59 

3.3.3 Summary 63 

3.4 Index Theorems and Spectral Flow 64 

3.5 Fermion Zero Modes in Other Theories 65 

3.6 A Model with Persistent Zero Modes 66 

3.7 Conclusions 67 

4 N = 1 Supersymmetric Abelian Cosmic Strings 70 

4.1 Introduction 70 

4.2 Supersymmetric Abelian Higgs Models and Spontaneous Symmetry Breaking 72 

4.3 Theory F: Vanishing Fayet-Iliopoulos Term 74 

4.4 Theory D: Non-vanishing Fayet-Iliopoulos Term 78 

4.5 The Super-Bogomolnyi Limit 79 

4.6 Concluding Remarks 80 

5 Cosmic Strings and Supersymmetry Breaking 82 

5.1 Introduction 82 

5.2 An SU(2) x U(l) Model 84 



iii 



5.3 SUSY Transformations and Zero Modes 85 

5.3.1 Abelian string 85 

5.3.2 Nonabelian string 86 

5.4 Soft SUSY Breaking 87 

5.5 SUSY breaking and an Abelian String Model 88 

5.6 Fermion Zero Modes after SUSY breaking 90 

5.6.1 U(l) Abelian models 91 

5.6.2 SU(2) x U(l) model 92 

5.7 Comments and Conclusions 92 

6 Massive Fermion Bound States 95 

6.1 Introduction 95 

6.2 Field Equations and a Positivity Condition 96 

6.3 Bound States in the Abelian String Model 98 

6.4 Bound States in Other Theories 99 

6.5 Summary 103 

7 Summary and Future Work 104 

7.1 Summary 104 

7.2 Future Work 106 

A An SO(10) Grand Unified Theory 110 

A.l Fermion and Gauge Fields 110 

A. 2 Higgs Fields 112 

A. 3 Fermion Masses 113 



iv 



List of Figures 

1.1 Finite temperature effective potential 5 

1.2 Plot of abelian string field profiles with m s = m v = 1 10 

2.1 SU{2) string field profiles, with rj G = A G = g = 1 and A' G = 5 23 

2.2 Plot of electroweak field profiles for abelian string 32 

2.3 Plot of electroweak field profiles for Wfj-string 35 

2.4 Variation of electroweak Higgs field at the centres of the abelian and W R 
strings 36 

3.1 The Dirac spectrum of cosmic strings with a zero mode and a very low lying 

bound state 64 

6.1 Spectrum of I = ... 3 fermion bound states 99 

6.2 Spectrum of I = 4 ... 7 fermion bound states 100 

6.3 Fermion bound state solutions in the abelian string model with mf/m s = 2.3. 101 



v 



Abstract 



In this thesis we investigate the microphysics of cosmic strings in non-minimal quantum 
field theories. In particular we consider theories in which fermion fields couple to the strings, 
and those with larger symmetry groups, such as grand unified and supersymmetric theories. 
By considering these extensions to the minimal model, we obtain a more realistic picture 
of the properties of cosmic strings. 

In considering grand unified theories, which have multiple phase transitions, we show 
that a cosmic string formed at one phase transition can cause the creation of another 
string-like solution at a later transition. This string-like solution will have many of the 
properties and implications of a normal cosmic string. We consider this effect for a general 
string solution, and illustrate it with a realistic 50(10) unified theory. As well as the 
usual abelian strings, this theory also contains more exotic string solutions. We consider 
both types of cosmic string. Separately, we examine the form of cosmic string solutions in 
supersymmetric theories, and the effect of soft supersymmetry breaking on them. 

We investigate the existence of conserved fermion currents in a variety of cosmic string 
models. We show that supersymmetry may be used to find the form of some solutions ana- 
lytically. We also derive an expression for the number and type of massless fermion currents 
in a general model. The existence of conserved currents can conflict with observations, so 
these results may be used to constrain models. We find the number of massless currents in 
the 50(10) and supersymmetric theories mentioned above. We show that currents present 
on a string can be destabilised by later phase transitions or supersymmetry breaking. This 
may allow any conflict that the current's existence has with observations to be avoided. We 
also examine massive fermion currents in a simple model, and determine the spectrum of 
such states. 
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Chapter 1 



Introduction 

1.1 Particle Cosmology and Cosmic Strings 

In the last twenty years there have been many significant developments in our understand- 
ing and approach to cosmology. In the past, the evolution of the universe was mainly 
studied through Einstein's equations of General Relativity. Modern cosmology also makes 
use of quantum field theory. This is particularly important when considering the early 
universe. According to the big bang model, the universe was once very small and at a very 
high temperature. Since classical approximations will certainly break down under these 
conditions, the investigation of quantum effects is vital. 

The combination of General Relativity and Particle Physics has led to a very accurate 
model known as the Standard Cosmology. This successfully predicts many phenomena, 
such as the expansion of the universe, light element abundances, and the uniformity of 
the cosmic microwave background. Despite its successes, there are still many unanswered 
questions. It is not clear how structures such as galaxies formed in a universe which seems 
to have been homogeneous and isotropic at very early times, nor how the universe became 
so uniform in the first place. There are also the questions of the origin of galactic magnetic 
fields, and why there is more matter than anti-matter in the universe. 

The Standard Model of particle physics does not provide answers to these and other 
questions. It is therefore necessary to go beyond the Standard Model and consider more 
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speculative ideas in quantum field theory, such as unification of the fundamental forces, 
or supersymmetry. Unfortunately, testing some of these ideas is beyond current particle 
accelerators. However, as well as predicting cosmology from quantum field theory, it is 
possible to use cosmological observations to make predictions about high energy quantum 
field theory. Study of the early universe can thus give insight into particle physics as well 
as cosmology. 

Topological defects (such as cosmic strings) are one example of quantum field theory's 
contribution to cosmology [1, 2]. They provide possible explanations for structure forma- 
tion [3], baryogenesis [4], cosmic microwave background anisotropics [5], and the origin of 
high energy cosmic rays [6, 7]. Once formed, a defect will not decay (unless it collides with 
another similar defect). Thus even though any defect formation will be restricted to the 
early universe, their effects will continue to the present. 

Topological defects form at phase transitions [8]. It is believed that the universe passed 
through several phase transitions shortly after it came into existence. These transitions 
reduced the symmetry of the field theory and gave mass to certain fundamental particles. 
This causes their interactions to be suppressed, and explains why they are not observed at 
everyday temperatures. Inside a defect the transition does not occur, so the laws of physics 
there resemble those of the universe before the phase transition at which they formed. 
This allows interactions to occur inside the defect which are heavily suppressed in today's 
universe. For example, if baryon violating processes were unsuppressed proton decay would 
occur more frequently [9] . This property of defects could provide an important window into 
the physics of the very early universe. 

In the past few years it has been realised that cosmic strings may have considerably 
richer microstructure than previously thought [10]. In particular, the presence of conserved 
currents in the spectrum of a cosmic string has profound implications for the cosmology 
of the defects. If the particles making up the current carry electric charge, and the string 
was formed at a mass scale of 10 16 GeV (a reasonable value for grand unification) the string 
could carry an electric current of up to 10 20 A [11]. There are many candidates for such a 
current, one class of which is introduced in section 1.6. 
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Another important consequence of cosmic strings is their gravitational effects. The 
mass per unit length of a string formed at an energy scale around 10 16 GeV would be 
10 22 g/cm. The gravitational field of the string is such that although matter is not attracted 
to a stationary string, it is attracted to the wake of a moving string. The evolution of a 
network of cosmic strings produced at high energy scales thus provides a possible origin 
for the seed density perturbations which became the large scale structure of the observed 
universe [3]. The gravitational effects of the string could also explain anisotropies in the 
cosmic microwave background [5]. 

The occurrence of phase transitions and defect formation is not restricted to quantum 
field theory. They also occur in condensed matter systems such as superconductors, 4 He 
and 3 He superfluids, and nematic liquid crystals [12]. Defects have actually been observed 
in these cases, and have been used to gain insight into the evolution of cosmological defects. 

1.2 Grand Unification and Symmetry Breaking 

In the Standard Model the weak and electromagnetic forces are unified at energies of order 
10 2 GeV. This idea can be extended to give a Grand Unified Theory (GUT) in which the 
strong force is also unified with the electroweak force. If this is the case, then it is believed 
that unification will occur at an energy scale of order 10 16 GeV. The unified theory will 
be based on a simple continuous Lie group, G. At high energies the vacuum state of the 
theory respects the full symmetry of the Lagrangian. Such a situation occurs in the early 
universe, when its temperature is extremely high. As it cools, the gauge theory undergoes 
a series of spontaneous symmetry breakings, until it becomes SU(3) C x U(1)q, which are 
the gauge groups of QCD and QED. This can be represented by 

G^H^ ► SU(3) C x SU(2) L x U(1) Y SU(3) C x U(1) Q , (1.1) 

where the first breaking occurs around 10 16 GeV. Unfortunately modern particle accelerators 
cannot hope to access energies above ITeV. The only place where suitable energies are likely 
to be reached is the early universe, so cosmology may provide evidence for unification. 
The spontaneous symmetry breaking mechanism works as follows: At high temperatures 
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the vacuum state, |0) GUT , is invariant under G, so g(x) |0} GUT = |0} GUT for all g(x) G G. 
At a lower temperature the theory's potential changes, and a different vacuum state, |0), 
becomes energetically favourable. This new state is only invariant under H, a subgroup of 
G, so the symmetry has been broken from G to H. 

Although the vacuum state is no longer invariant under G, the theory itself must still be. 
Thus if |0) is a minimal energy vacuum state, so is g |0), hence the vacuum is degenerate. 
The vacuum manifold then consists of all the distinct states of the form g | 0) . Since g\ |0) = 
92 |0) if 9i 1 g2 G H, the vacuum manifold M. is equal to the coset space G/H. 

The different vacuum states can be labelled by the expectation value of a scalar Higgs 
field. In a general GUT there will be several Higgs fields, each associated with a subsequent 
phase transition at which symmetry is broken. As a simple example consider a model with 
G = U(l) and H = I. It has the Lagrangian 



where <\> is the (complex) Higgs field, D^cf) = (d^ — ieA^) 4> and F^ u = d^A u — d u A^. V 
is the potential energy of <p. It is equal to (A/4)(|</>| 2 — rj 2 ) 2 at zero temperature. When 
thermal corrections are added, its minimum varies with temperature. Above some critical 
temperature T c the vacuum energy is minimised by 4> = 0. Below this temperature the 
minimum has <j) ^ 0. T c is the temperature at which the phase transition occurs. If T is 
the temperature of the universe, a suitable V is 



Thus the vacuum states have = for T > T c and <f> = r)e ia y/l — T' 2 /T^ = constant for 
T < T c (see figure 1.1). <fi = r]e ta when T = 0. The above potential gives a second order 
phase transition. 

The gauge and fermion fields of a GUT will couple to the Higgs fields. When a Higgs 
field gains a non-zero vacuum expectation value (VEV) , it gives masses to some of the gauge 
and fermion fields via this coupling. The masses from a GUT scale phase transition will be 
enormous, and so the corresponding particles and their interactions will not be observed. 
This explains why we do not observe significant proton decay in today's universe. However, 



C = (D tl 4>)*(D tJ -(f>) - -F^ - V{4>) , 



(1.2) 




(1.3) 
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Figure 1.1: Finite temperature effective potential. 

such forces must have been unsuppressed at some time in the early universe, since the uni- 
verse could not have reached its current state without them. Experimental measurements 
of the variation of gauge couplings suggest that they converge at high energies [13]. This 
is further circumstantial evidence for unification. 

In a grand unified theory there will be one unified force, and one unified fermion field. 
The idea of supersymmetry (SUSY) is a natural extension of this, in which the gauge and 
Higgs fields are unified with the fermion fields [14]. If the strengths of the fundamental 
forces are carefully extrapolated to high energies, they do not quite converge, as would be 
required in a GUT. This problem is solved if SUSY is added to the theory, suggesting it 
does have physical relevance. Furthermore, supersymmetry and grand unification form a 
central part of superstring theory, which is the most promising theory to describe gravity 
on the quantum scale. Since any full theory of everything includes quantum gravity, it 
seems likely that SUSY must play some role in a realistic quantum field theory. Another 
motivation for supersymmetry is its solution to the hierarchy problem [15]. We expect the 
high energy scale GUT fields to couple to the low energy scale Higgs field of the Standard 
Model. The mass of this Higgs field receives quadratically divergent radiative corrections 
from these couplings, which will generally lead to conflict with experiment. The symmetry 
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of a supersymmetric theory ensures that the contributions from the fermion and boson 
fields cancel each other exactly, avoiding the problem. 

1.3 Topological Defects 

Phase transitions would have occured in the early universe as it cooled down. Because 
information can travel no faster than the speed of light, (j) would gain different vacuum 
expectation values in parts of the universe that were not in causal contact [8]. It is ener- 
getically favourable for these variations to disappear. However if M has any non-trivial 
homotopy groups it is possible that topologically stable configurations with (f> not con- 
stant will form. These are called topological defects. The type of defect depends on the 
non-trivial homotopy group. 

If tto(M) / / then M is disconnected. Suppose that <p 6 Ma at x = od (cartesian 
coordinates) and <p G Mb at x = — oo, where Ma and Mb are disconnected components 
of M. By continuity there must be a region where <p M. This region is called a domain 
wall. In three dimensions it will be a two dimensional surface. 

In the model given by (1.2), iri(G/H) = Z / /, so M contains non-contractible loops. 
In two dimensions, a configuration with 4> = rje md at infinity (with n any integer) cannot be 
continuously deformed to 4> = constant, and so is topologically stable. Regularity will force 
4> to leave M in some region (see next section). Such defects are called cosmic strings, as 
in three dimensions they are one dimensional. For topological reasons it is not possible for 
cosmic strings to have ends, so they must either be infinitely long, or a closed loop. Suppose 
a cosmic string did have an end. Consider the variation of the Higgs field along a closed 
path around a string. This variation corresponds to a non-contractible loop in M. By 
sliding it off the string we can continously change this path to one which corresponds to a 
contractible loop in M. This is equivalent to continously moving between the disconnected 
parts of H, which is not possible. Hence strings cannot have ends. 

In three dimensions, if iT2(G/H) / / the corresponding point-like defects are called 
monopoles. In this case there are mappings of S 2 to itself which cannot be continously 
deformed to the identity. Thus if M is equal to S 2 , and (p at r = oo is such a map, the 
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resulting configuration will be topologically stable. As with strings, there will be a region 
in which <fi M.. 

The size of a defect can be estimated by balancing the potential and kinetic terms of 
its Lagrangian (1.2). The potential energy is of order Ar/ 4 inside the defect. If 5 is the 
width of the defect, the kinetic terms (d(j)) 2 are of order (f]/5) 2 . Equating these terms gives 
5 ~ l/(V\rj) ~ m^ 1 , where m s is the mass of the Higgs field. 

The cosmological consequences of domain walls and monopoles strongly conflict with 
observations. Domain walls will come to dominate the energy density of the universe [16], 
and the predicted monopole density is unacceptablely high [17]. Thus any GUT which 
predicts them must also have a mechanism which ensures most of them are destroyed. 
The properties of cosmic strings do not conflict with observations (see section 1.7). Fur- 
thermore they provide explanations of many phenomena, so cosmic strings are the most 
cosmologically significant defect. 

1.4 The Abelian Cosmic String Model 

As it stands, the string solution described in the previous section has a non-vanishing 
covariant derivative which gives an infinite contribution to the energy. This is avoided by 
having a non-zero gauge term. The resulting solution is a cosmic string. More precisely 
(taking T = 0), it has the form 



In order for the solution to be regular at the origin /(0) = a(0) = 0. If the solution is to 
have finite energy, f(r) and a(r) must tend to 1 as r — ► oo. Substituting (1.4) and (1.5) 
into the theory's field equations gives 



4> = Vf(r)e 



(1.4) 



A 




(1.5) 




/ = ^ 2 (/ 2 " 1)/ , 



(1.6) 



a" - - = -2eV (1 - a)/ 2 . 



(1.7) 
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The resulting string is the well known Nielsen-Olesen vortex [18]. It turns out that / and a 
take their asymptotic values everywhere outside of a small region around the string. Thus 
\<f>\ is constant and is pure gauge away from the string. The sizes of the regions in which 
the magnetic field is non-zero and \4>\ is not constant (r v and r s ) are roughly the inverses 
of the masses of the corresponding particles. Thus if m s and m v are the Higgs and gauge 
field masses, 

r^ 1 «m s = V~\r) , r^ 1 m m v = \pier] . (1.8) 

By rescaling r, the e and rj dependence can be removed from (1.6) and (1.7). The form of 
the solutions then only depends on only one parameter: = (m s /m v ) 2 . 

The above solution is a point-like defect in two dimensions. If it is extended to three 
dimensions it will take the form of an infinite line, or a closed loop. The defect is now one 
dimensional and string shaped, hence its name. There is a non-zero magnetic field inside 
the string. Its flux is 

<S> B = / A^ = . (1.9) 

J r=oo & 

Many of the observational consequences of cosmic strings arise from their gravitational 
effects. Defining fi to be the mass per unit length of the string, the size of these effects is 
proportional to the dimensionless quantity G\i ~ (T c /m p i) 2 , where G is Newton's constant 
and m p i is the Planck mass. We find that \i ~ i] 2 , so for a grand unification scale string 
G/j, ~ 1CT 6 . This is the right magnitude to explain the fluctuations in the cosmic microwave 
background and the matter distribution of the observed universe. 

The spacetime around a string is approximately conical. It resembles a flat space with 
a wedge of angular size 8irGfi [19] removed, and the two faces of the wedge identified. If 
two objects travelling on parallel paths pass each side of a string, they will begin moving 
towards each other despite having experienced no force. The result of this effect is that as 
a string moves about, it increases the density of the regions through which it passes. This 
is an example of how cosmic strings provide a mechanism for structure formation. 

At the centre of an abelian cosmic string (f> is zero. Since the breaking of G to H is 
caused by <j> being non-zero, this means that G is not broken inside the string. The result is 
also true for other types of defects (monopoles and domain walls). Thus the grand unified 
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theory is restored inside the string, even though the universe has cooled to temperatures 
that would usually break it. 

The above ansatz (1.4,1.5) describes an infinite number of string solutions since the 
winding number, n, can be any non-zero integer. Although they are all topologically stable 
with respect to the vacuum and other strings, it is possible for |n| > 1 strings to decay by 
splitting into several strings with smaller winding numbers [1]. 

If (3 > 1 (where (3 = (m s /m v ) 2 ) this does happen, since the force between the strings 
is repulsive. The reverse is true when (3 < 1. In this case strings will tend to combine to 
produce a single string with a larger winding number. The situation is similar in supercon- 
ductors. In this case (3 < 1 and [3 > 1 correspond to type I and type II superconductors 
respectively. 

When (3 = 1 there is no force between strings. It is then possible to reduce (1.6) and 
(1.7) to first order equations: 

f = \n\ f -{l-a) , (1.10) 

H^ = ^ 2 (i-/ 2 ). (i.ii) 
1.5 Determination of the String Field Profiles 

Although full analytic solutions of the string field equations have not been found, their 
asymptotic forms are easily determined. For large r, f(r) ~ 1 and (1.7) can be reduced to 
the equation for a modified Bessel function. The solution is 

l-a(r) ~rA"i(m v r) . (1.12) 

When (3 > 4 the gauge term in (1.6) determines the behaviour of f(r) [20]. Otherwise it 
can be dropped and 

1 - f(r) ~ K (m s r) . (1.13) 
The forms of f(r) and a(r) at small r are 

f{r) ~r>l , a{r) ~ r 2 . (1.14) 
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2 4 r 6 8 10 

Figure 1.2: Plot of abelian string field profiles with m s = m v = 1. 

It is possible to find solutions of the string field equations numerically. One suitable 
method is shooting. In this method, values for the free parameters in the small r solution 
(1.14) are guessed, and the equations are then integrated out to large r. The values of f(r) 
and a(r) are then compared with the required values. This is repeated several times, with 
the choice of free parameters being adjusted each time to get closer to the required solution 
at large r. 

Rather that trying to match boundary conditions at infinity, a large finite value of r 
will do. In fact, as well as the required solution, (1.6) and (1.7) have an exponentially 
increasing large r solution. Thus the accuracy of the computer will severely limit the range 
of r which can be sensibly used. Alternatively it is possible to shoot from both small and 
large r, using (1.12,1.13), and require that the solutions meet at some point in the middle. 

A better, but more complex and memory intensive approach is relaxation. The idea 
is to start with an approximate solution for all r, and then see how closely it satisfies the 
field equations. The difference from the required values of the field equations is then used 
to calculate an adjustment to the guessed solution. This is repeated several times until 
the modified guess satisfies the field equations closely enough. Unlike shooting, the whole 
solution is stored in the computer's memory, so this approach uses more resources. In the 
case of cosmic string field equations, it is very reliable. An example of solutions obtained 
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by relaxation is shown in figure 1.2. 



1.6 String Superconductivity 

In the core of a cosmic string the Higgs and gauge fields do not take their usual values. 
It is possible that other fields will also have different VEVs, resulting in the string having 
additional properties. The first example of this was discussed by Witten [11]. He considered 
an abelian string theory with an extra scalar field, a, whose usual VEV is zero. For a suitable 
choice of potential this scalar field gains a non-zero VEV in the core of a string (due to the 
variation of 4>). Since a is charged, its non-zero VEV will break electromagnetism. If z and 
t dependence are added to a, the resulting solution is a superconducting current, which is 
conserved. Since a is electrically charged, the current will have long range effects. 

It is also possible to have currents made up of fermions or gauge bosons. Gauge boson su- 
perconductivity only occurs in more complex nonabelian theories, and will not be discussed 
here. Fermion superconductivity can occur in simpler models, such as the abelian model 
discussed in section 1.4. The fermion superconductivity arises from Yukawa couplings, and 
unlike scalar boson superconductivity, it is not necessary to use a specific potential, so 
fermion superconductivity is more generic. Strictly, cosmic strings with fermion currents 
are not superconductors since electromagnetism is not broken inside the string. The strings 
act as perfect conductors, since the current flows without resistance. However, the term 
'superconductivity' is often used for both cases. 

If charged currents can exist on a cosmic string, they will be generated when the string 
passes through an electric or magnetic field. Currents can also be generated by interaction 
with the plasma, particularly when the string forms, or by string collisions. It has been 
suggested that internal phase transitions occuring on the string could generate currents 
too [21]. The presence of a current will alter the evolution of a string network. Decaying 
currents may explain observed high energy cosmic rays [7, 11], and could also provide a 
mechanism for baryogenesis [22]. Charged currents on a fast moving string can create 
Shockwaves in the plasma, which has implications for structure formation [23]. It is also 
possible that currents could stabilise loops of cosmic string [24]. Such configurations are 
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called 'vortons'. If they form at high energy scales they can have dramatic consequences 
(see next section). Vortons formed at low energies may provide a dark matter candidate [1]. 

The maximum fermion current on a string is restricted by the mass of the corresponding 
particle off the string. If the current's momentum exceeds this it can escape from the string. 
The ease with which the current then escapes will depend on the curvature of the cosmic 
string. For a string with radius of curvature R, the following bound on the current is 
obtained [25] 

Jmax ~ em 2 R . (1-15) 

In fact the maximum value the current will reach is likely to be less than this, and will 
depend on the details of the model. 

Consider an extension of (1.2) to include a two-component fermion, with charge 1/2. 
The extra terms in the Lagrangian will then be 

Acrmions = ^lO^D^ - - [igy^V + (h- C.)] , (1.16) 



where cr M = (—1, a 1 ), D^ip = [d^ — ^ieA^j ip, and ip c = ia 2 (^) T is the charge conjugate of 
■ip. This gives the field equations 



V 



dr + ide + n^] d z + d t 



ip - m f /(r)e in V = , (1.17) 



d z -d t d r -lde-n^ 



2r 

where the expressions (1.4,1.5) have been substituted for <p and A^, and mf = gyV- 

Non-trivial solutions of (1.17) with only r and 9 dependence exist. It has been shown 
that there are |n| such solutions which are normalisable, in the sense that j\^\ 2 d 2 x is 
finite [26]. They have zero energy, and are referred to as zero modes. If n = 1 the single 
solution can be found analytically, 

V(r, 9) = Q exp (- £ m { f(s) + ^ds^j . (1.18) 

For higher \n\ it is not generally possible to determine all such solutions analytically, al- 
though some of their properties can be found. Like (1.18), they are all eigenstates of a 3 . 
Their eigenvalues are +1 if n > 0, and — 1 if n < 0. All the solutions decay exponentially 



12 



outside the string, and so are confined to it. They can be regarded as fermions trapped on 
the string. 

The solutions can be extended to include z and t dependence. This is achieved by 
multiplying ip by a(z,t), which satisfies (d z =p dt)a = 0, depending on whether a 3 ip = ±"0- 
Thus the trapped fermions move at the speed of light, in the ±z direction. Other fermions 
which couple to <ff can be added to the theory. </>* has the opposite winding number to cj>, 
and the fermion currents flow in the opposite direction. Considering both sets of fermions, 
currents can flow in both directions. 

Although ip(r,9) was easy to find for n = 1, this is not generally the case for other 
winding numbers, or for more complex theories with several fermion fields (but see chap- 
ter 4). More realistic grand unified theories will have several Higgs fields, and possibly more 
complex string solutions than the abelian case discussed in section 1.4. It is still possible 
to determine the existence of fermion zero modes in such cases. This is done for a general 
theory in chapter 3. A specific grand unified theory is also considered there. 

In supersymmetry the fermion and boson fields are related by a symmetry. This means 
that the properties of a cosmic string (which is a bosonic object) can be used to make 
predictions about the fermion solutions on the string. This is a useful idea, and it is 
investigated in chapters 4 and 5. 

1.7 Evolution of Cosmic String Networks 

Although cosmic strings form in the early universe, many of their observable effects will 
take place in the later universe. It is therefore important to know how a network of cosmic 
strings will evolve. 

When considering a network of cosmic string it is useful to define a characteristic length 
scale £. There are various possible choices of £, for instance £ 3 can be the volume of space 
which contains cosmic string of average length £. Initially £ will grow faster than the size 
of the universe [2]. This cannot continue for long since causality implies that £ < t. 

Logically, £ can do one of two things after this point. It can either approach a scaling 
solution where £/t is constant, or it can grow less quickly, so £/t decreases. If £ grows 
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more slowly, the strings will come to dominate the energy density of the universe. This 
is sometimes referred to as 'overdosing' the universe. This would substantially alter its 
evolution. Such evolution is strongly ruled out by observations. 

The most obvious way for a string network to evolve is to just stretch with the expansion 
of the universe. Unfortunately this implies £ oc y/i, and so the strings dominate the energy 
density. To avoid this, there must be some mechanism for transferring energy away from 
the string network. 

Strings can lose energy by radiating particles. Unfortunately strings are not usually 
charged, and their interactions with other fields are likely to be weak. Strings formed at 
high energy scales will emit significant amounts of gravitational radiation, although not 
enough to solve the problem with the network evolution. 

For topological reasons, strings are either infinite or closed loops. String loops provide 
a solution to the energy density problem. When a string intersects itself, a loop will 
break off. The loop will then start losing energy by gravitational radiation. It will begin 
to contract, and will eventually disappear. Unless this loop rejoins the string network, its 
contribution to the network energy density is lost. Smaller loops are unlikely to do this, and 
loop formation provides a substantial enough energy loss mechanism to ensure the string 
network has a scaling solution. The decay of string loops can also explain the observed 
baryon asymmetry [27]. 

If conserved currents exist on the strings, such as those discussed in section 1.6, the above 
solution fails. It is possible for string loops to be stabilised by the angular momentum of 
the trapped charge carriers [24]. Such stable loops are called vortons. The loops do not 
decay, so they will continue to contribute to the energy density of the cosmic string network. 
Since it no longer has a sufficiently strong mechanism of energy loss, the string network will 
no longer be able to reach a scaling solution. The universe will then become dominated 
by vortons. The possibility of this happening allows the underlying particle physics to be 
cosmologically constrained [28] . The study of fermion (and other) currents will give insight 
into the possibility of vorton formation. 

When strings collide they will intercommute. It is because of this that a loop is formed 
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when a string intersects itself. Some strings, such as type I (rrig/m^. < 1) abelian strings, 
do not intercommute [29]. In order for such strings to be cosmologically viable, they must 
either form late in the universe (and so not have time to dominate it), or possess a new 
mechanism for energy loss. Because of this, type II strings are likely to be more physically 
significant than type I strings. 

Similar arguments can be applied to domain walls and monopoles. There is no equivalent 
of loop formation for domain walls. These defects would be sure to dominate the energy 
density of the universe, hence they are ruled out [16]. The monopole density evolves the 
same way as the matter density. However, monopoles have long range effects which strongly 
conflict with observations [17]. Thus they are ruled out too. 

1.8 Overview 

Although a simple U(l) model is useful for illustrating the existence of strings and their 
basic properties, a more realistic theory is needed to get accurate phenomenological results. 
Furthermore, strings in such theories will have a richer microstructure, which can give rise 
to additional properties and phenomena. It may also alter, or even remove, some of the 
properties suggested by the simplest models. In this thesis we will consider cosmic strings 
in grand unified and supersymmetric theories. Perhaps the most natural place for cosmic 
strings to arise is in unified theories, so it is important to consider the implications of their 
more complex symmetry breaking for strings. Many credible extensions of the Standard 
Model involve supersymmetry, and so its effects should also be examined. 

In chapter 2 a grand unified theory which has cosmic strings is outlined. Its unifying 
symmetry group is SO( 10). There is a wider range of string solutions in this theory than 
the U(l) model. Many of these have a more complex internal structure. Instead of just one 
winding number, their Higgs fields will have several. We show that a cosmic string can affect 
subsequent phase transitions, forcing other Higgs fields to take string-like solutions. We 
investigate this effect in the 50(10) theory. Generalisation of the results to other related 
theories is also discussed. 

As was mentioned in section 1.6, precise determination of fermion zero modes is rarely 
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possible. Their existence can be determined by examining the field equations, but this is 
time-consuming for more complex theories. In chapter 3 an index theorem is found which 
gives the number of zero modes for a general theory. Whilst index theorems have been 
found previously, they only give the difference between the number of left and right moving 
currents. Ours is more general, and gives the total number of massless currents. By apply- 
ing the theorem before and after a phase transition we can investigate the fate of fermion 
zero modes during sequences of phase transitions in a variety of models. Depending on the 
couplings that the breaking introduces, the zero modes may be destroyed and the supercon- 
ductivity of the string removed. Vortons will then dissipate, relaxing the constraints on the 
theory. We discuss the features of the theory that are required to produce this behaviour 
and consider the implications of spectral flow. We apply the theorem to the unified theory 
of chapter 2. It is applied to other theories in sections 4.3 and 5.3. 

In chapter 4, the microphysics of supersymmetric cosmic strings is discussed. For sim- 
plicity two N = 1 supersymmetric abelian Higgs models are considered. The vortex solu- 
tions are found, and it is shown that the two simplest supersymmetric cosmic string models 
admit fermionic conductivity. In a SUSY theory, fermion and boson fields can be trans- 
formed into each other. This allows string solutions with non-zero fermion fields to be found 
analytically. These solutions are fermion zero modes, and are found to first order explicitly. 
We note that this constrains all supersymmetric grand unified theories with abelian strings. 

In chapter 5 we extend the results of the previous chapter to a model with abelian 
and nonabelian strings. The strings in this model have some resemblance to those of the 
50(10) model discussed in chapter 2, and so may give some insight into the properties 
of SUSY GUT strings. We give the string solutions, and find analytic fermion zero mode 
solutions using SUSY transformations. We consider the effects of soft supersymmetry 
breaking on these cosmic string solutions, as well as those in the previous chapter. We also 
examine the implications of SUSY breaking for the fermion zero modes. Soft SUSY breaking 
terms are those which break SUSY without giving quadratically divergent corrections to 
the electroweak Higgs field. 

So far we have only considered massless fermion currents, arising from zero modes. In 
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chapter 6 we look for massive currents arising from bound states. We show that there 
are no space-like fermion currents in any model. In contrast to the null (or light-like) 
currents, it is difficult to determine the existence of the time-like currents analytically, so 
we use numerical methods instead. We determine the spectrum of fermion bound states 
and currents in the abelian string model. We also speculate about similar states in other 
theories. 

Finally, in chapter 7, the various results are summarised, and possible future work is 
discussed. 
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Chapter 2 

Microphysics of SO(10) Cosmic 
Strings 

2.1 Introduction 

As was discussed in chapter 1, many significant cosmological properties of cosmic strings 
arise from their microstructure. In particular, at subsequent phase transitions the core of 
the cosmic string acquires additional features. For example, the string can cause electroweak 
symmetry restoration in a much larger region around it, proportional to the electroweak 
scale itself [10, 30]. This microphysical structure has been used to provide a new scenario for 
electroweak baryogenesis [31], and to investigate the current-carrying properties of cosmic 
strings [32, 33]. 

Previous work considered the simplest extension to the Standard Model that would 
allow the formation of strings. A U(l) symmetry, whose breaking produced an abelian 
string, was added to the usual Standard Model symmetries. The resulting theory had two 
coupling constants, of arbitrary ratio. It was shown that if the ratio was large enough, the 
electroweak Higgs field would not only be zero at its centre, but would also wind like a 
string. Whether this is likely to happen with phenomenological strings can be found by 
considering a realistic grand unified theory, where there is less arbitrariness. 

By using a larger gauge group it is also possible to consider the effects of nonabelian 
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strings, which could not occur in the theories considered in refs. [10, 30]. Nonabelian 
strings have significantly different behaviour to abelian strings, since the associated string 
generators do not all commute with the Standard Model fields, or the other gauge fields. 
It is thus necessary to approach them in a slightly different way. 

In this chapter we examine these issues in detail for strings formed in a realistic grand 
unified theory (GUT) based on 50(10). In section 2.2 the theory to be used is outlined. 
The possible strings that form in it prior to the electroweak phase transition are discussed 
in section 2.3. The effects of a cosmic string on the subsequent phase transitions of a 
general theory are considered in section 2.4. In section 2.5 we consider these effects in 
greater detail for the electroweak phase transition of the 50(10) GUT, and examine the 
various electroweak string solutions in section 2.6. The form of the electroweak fields, and 
numerical solutions of the corresponding field equations, are found. In section 2.7 some 
other, simpler symmetry restorations occuring in the theory are discussed, in particular 
that of the intermediate SU(5) symmetry. Although one specific theory is considered, 
many of the results generalise to other theories. The implications of our results for such 
theories are discussed in section 2.8. In section 2.9 we summarise our results and discuss 
the conclusions. 



2.2 An SO(10) Grand Unified Theory 

A realistic GUT which has a symmetry breaking pattern which produces strings is 5O(10). 
Its properties have a reasonable agreement with physical results. Consider the symmetry 
breaking 

5O(10) SU(S) x Z 2 

^ SU{3) C x SU(2) L x U(\) Y x Z 2 

^ SU(3) C x U{l) Q x Z 2 , (2.1) 

where 3>n transforms under the N representation of 5O(10). The actual grand unified 
gauge group is Spin(W), the covering group of 5O(10), but for simplicity the symmetry 
breaking is shown in terms of the Lie algebras. The discrete Z 2 symmetry left by the <I>i26 
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Higgs field leads to the formation of a variety of cosmic strings. Comparison of the effects 
of the various symmetry breakings is simplified by expressing everything in terms of the 
same representation. Since 126 + 10 = (16 x 16)5 this is possible. Conveniently, 16 is 
also the representation that acts on the fermions. The 50(10) fermions consist of the usual 
standard model fermions, plus a right handed neutrino. The fermionic part of the theory 
is then expressed in terms of the left-handed fermions and the charge conjugates of the 
right-handed fermions. The fields of the theory are discussed in more detail in appendix A. 

The maximal subgroup of 50(10) is actually SU(5) x U(l)p, and P can be used to 
decompose 50(10) into representations of 5C/(5) 

16 ^ 1 5 + 10i + 5_ 3 , (2.2) 

where the subscripts are the eigenvalues of P. 126 and 10 can be similarly decomposed by 
considering symmetric products of 16. 

126— ►I10 + ... , 10— >5_ 2 + 5 2 . (2.3) 

P can also be used to describe the non-trivial element of the discrete symmetry of (2.1), 
which is d = exp(27riP/10). 

Defining to be the usual constant vacuum expectation value of <3>n, 3>i26 nas a 
magnitude of r] G , which is of order 10 16 GeV. It is in the lio component of 126, and so 
must be equal to rj G (ei x ei), where ei is in the I5 component of the 16 representation 
(the corresponding field in the fermion representation is the charge conjugate of the right- 
handed neutrino). ^Jjf is made up of both the chargeless components of 10. If H® and 
H® are the chargeless components of 5_2 and 5 2 respectively, ^iq = VdH^ + r] u H®. Since 
10 is contained in (16 x 16)5, $10 can be expressed as a sum of symmetric products of 
components of 16s. The 45 is contained in 16 x 16. $45 's effect on the form of the cosmic 
string solutions is far less significant than the other Higgs fields, so we will ignore it for 
now. 

Strings can form at the first 5O(10) — ► SU(5) x Z 2 symmetry breaking. In this case 
$126 is not constant, and takes the form e ldTs Q^si 7 ') ■ ^126 ^ s independent of 9, and satisfies 
the boundary condition $126 (°°) = ^126 ■ * s ma de up of the broken generators of 5O(10), 
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and must give a single-valued <l?i26- If e 2mTs = U x d (for some U in SU(5)), then while 
<I>126 will be single valued, it will not be topologically equivalent to $126 = constant, and 
so the string will be topologically stable. If e 2mTs = U x I the string is not topologically 
stable, but may have a very long lifetime, and so still be physically significant [11]. 
The Lagrangian of the system is 

C = {D^ 12& )%D^ 126 ) + {D^ w )\D^ w ) + {D^ Ab )\D^ b ) 

~ \F% F ^ a - ^(*126, $45, $10) + Acrmions , (2.4) 

where = - ^igA^ and F^ u = d^A u - dyA^ - \ig\A^A v \. There are 45 gauge 
fields in all, most of which acquire superheavy masses and so are not observed at everyday 
temperatures. They consist of the usual standard model fields, with the VF-bosons denoted 
by Wl; Wft, which are right-handed versions of the W£, coupling right handed neutrinos 
to electrons; some leptoquark bosons: , X^ , X^, where Yf 1 and Xf 1 are SU(5) gauge 
fields; some more general gauge fields X'^ and , which couple quarks to leptons and 
different coloured quarks; and a fifth uncharged field, B' . The index i takes the values 1, 
2, 3, and is related to colour. Two uncharged diagonal fields, Z' and B, are made up of 
orthogonal linear combinations of and B'. Linear combinations of B and produce 
the Z boson and the photon, A. At the first symmetry breaking Z', W^, X'^, Y(^, and 
Xg f are all given superheavy masses. The second stage gives high masses to Xf, Y^, and 
additional masses to X-^, Y^ and X^. Finally the third symmetry breaking gives 
masses to W^and Z, with further masses being given to the Z', W^, Y^ and Y(^ fields. 



2.3 GUT Strings 

Neglecting fermions, the Euler-Lagrange equations obtained from (2.4) are 

D,D^ = -|^, (2.5) 

i 

(D.F^r = -glmX;^)*^) , (2.6) 

i 

where D tl F tlv = d^F^ - if [A^, F^]. At high temperatures V is such that $126 is the only 

non-zero Higgs field. (2.5) and (2.6) have various cosmic string solutions. The different 
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solutions correspond to different choices of T s . In the 50(10) — ► SU(5) x Z2 symmetry 
breaking, 21 of 50(10) 's 45 generators are broken, and T s will be a linear combination of 
them. One of them, P, corresponds to the U{\) symmetry not embedded in SU(5). The 
corresponding string is abelian, and has the solution 

$126 = f{r)e m0T ^Xf, , 
2a(r) 

Aq = n — —T s , An = otherwise , (2.7) 

gr 

where T s , the string generator, equals P/W, and n is an integer. The non-zero gauge 
field is required to give a zero covariant derivative, and hence zero energy, at r = 00. 
It corresponds to a non-zero Z' field. (2.7) can be simplified using Pe\ = 5ei, to give 
e m0T s$ yac = e m0 $ yac Substituting (2.7) into (2.5) and (2.6) gives the Nielsen-Olesen 
vortex equations, as would be expected. Regularity at the centre of the string, and finite 
energy due to a vanishing covariant derivative and potential at infinity, imply the boundary 
conditions /(0) = a(0) = and /(oo) = a(oo) = 1. 

The situation for the other generators is more complicated. For a general string genera- 
tor T s , the left and right hand sides of (2.5) are proportional to T^^e an d $126 respectively, 
which in general are not proportional [34]. Thus the solution (2.7) will not work. This is 
resolved by expressing $||6 i n terms of the eigenstates of T s 2 , to give = Em 4>m, where 
T 2 4> m = m 2 4> m . Since T s is Hermitian, m 2 will be positive and real. A suitable string 
solution can now be constructed 

$126 = e meT ^fm(r)^m, 
m 

2a(r) 

A e = n — —T s , A u = otherwise . (2.8) 

gr 

In order for $126 to be single valued the various m must all be integers, and T s <fio must 
be zero. The boundary conditions on a and f m will be the same as those for (2.7), except 
that fo need not be zero at r = 0. The simplest examples of such solutions occur when 
T s 2 ei = |ei, in which case 

4>i = y [(ei x ei) + 4(T sei x T s ei)] , (2.9) 
4>o = y [(ei x ei) - 4(T sei x T^)} . (2.10) 
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Figure 2.1: SU(2) string field profiles, with rj G = A G = g = 1 and A' G = 5. 

T s ei and ei are orthogonal, so 0o and 0i are orthogonal. In this case only part of the 
Higgs field winds around the string. This type of string was first suggested by Aryal and 
Everett [35], and has been examined in detail by Ma [34]. It turns out to have lower (about 
half as much) energy than the abelian string (2.7). This is because the Higgs field is not 
forced to be zero at the string's centre, which reduces the contribution to the energy from 
the potential terms. Also since $126 varies less, the covariant derivative terms are smaller. 

Of course, such vortex-like solutions are only topological strings if e 2m 8 is not contained 
in SU(5). If n is even this is not the case, and the solution is topologically equivalent to 
the vacuum. Similarly, odd values of n are all topologically equivalent to each other, so 
there is only one topologically distinct type of string of this form. Strings with higher n 
can unwind into strings with lower n. The same is true of the abelian string. However it 
is possible that the lifetime of an n > 1 string will be very long, so in a general theory all 
values of n should be considered. 

As shown in [34] the most general potential reduces to a different form to that of the 
abelian case, and leads to these equations for a and the / m 's 



fi + 7 - n 2{ ±^h = f [Ac (fl + fo - 2) + A' G (ft - /o 2 )] h , (2.12) 




Ac (/i 2 + /o 2 -2 



:)-A G (.A 2 -/o 2 )].f( 



(2.11) 
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a"-^ = -g 2 V 2 G (l-a)f 2 1 . (2.13) 

A G and A' G are such that /i(oo) and /o(°°) will both be 1. A numerical solution of the 
above equations is shown in figure 2.1. Since /o / at the string's centre, symmetry is still 
broken there (although the VEV of <3?i26 is lower). The value of /o(0) depends on the ratio 
of A G and A' G . If A' G > A G , / (0) < 1 while A' G < A G gives / (0) > 1. If A' G = A G , f = 1 
everywhere, and (2.12,2.13) reduce to abelian string field equations. 
The corresponding equations for the abelian string (T s = P/W) are 

/" + £ - n 2 ^^/ = ^ G A G (/ 2 - 1)/ , (2.14) 

a"-^ = -5gW G (l-a)f 2 . (2.15) 

The above nonabelian strings are in fact all SU(2) strings. There are other more com- 
plicated possibilities, for which T s 2 ei is not proportional to ei, but none of these are 
topologically stable. We shall only consider topologically stable strings, and the closely 
related solutions with higher winding numbers. We will consider the strings corresponding 
to each of the broken generators. These are all equivalent under SU(5), but not under 
SU(3) C x U(1)q, so they will be distinct after the electroweak symmetry breaking. Apart 
from the abelian string the four cases correspond to non-zero Wg, X^, X'± and Yf* fields. 
Under SU(3) C x U(1)q any linear combination of generators can be gauge transformed 
into any other combination, thus they are equivalent. The same is true for the other gen- 
erators, so there are just 5 distinct types of string at low temperatures. The 4 nonabelian 
strings can be labelled by their gauge fields. Under U(1)q, nonabelian strings with winding 
number — n are gauge equivalent to ones with winding number n (for any choice of T s ), so 
it is sufficient to consider only n > strings. 

Since all the nonabelian strings are gauge equivalent under SU(5), they have the same 
energy, which is about half that of the abelian string [34]. Later phase transitions will alter 
the string solutions and remove this degeneracy. 
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2.4 Cosmic Strings and Multiple Higgs Fields 

Before examining the effect of the other phase transitions on the 50(10) cosmic strings, we 
will consider a general GUT. Suppose it has the symmetry breaking 

G Gi x D ^ G 2 x D ■ ■ ■ G n _i xfl^G„xD (2.16) 

where D is a discrete group, such as Z^. Define 3>J ac to be the usual constant vacuum 
expectation value of <3?j, and = — i^A^ to be the covariant derivative. Because of the 
discrete group cosmic strings can form at the first phase transition. Away from the string 
core the solution will take the form 

$o = Ro(eW ac , A = —R o \e)d e Ro(6) , (2.17) 

igv 

where Rq(8) G Gq. We will choose Ro(0) = I. Since <E>o is single valued, R must obey 
i?o(27r)^ ac = $|f c , thus i? (2vr) = Ui x d, where U\ G G\ and d G D. If the string is to be 
topologically stable d must not be equal to /. The gauge field ensures a vanishing covariant 
derivative, and thus finite energy. 

At the next phase transition $i gains a VEV. If its covariant derivative is also to vanish, 
<J>i must wind like <!>o, so 

$i = i? (6W ac • (2-18) 
Unfortunately if Ro(2-!r)^] ac / ^Y^i then $i is not single valued. The same is true of gauge 
and fermion fields who get their masses from $i. Therefore this ansatz is unacceptable. 

If the product G\ x D is a direct product, it is possible to alter the above ansatz to give 
a single valued <E>i, without changing $ - This is achieved by applying an extra (singular) 
gauge transformation that has no effect on 3>o> but ensures the values of $1 at and 2n 
match up. A corresponding gauge term is then added so the covariant derivative vanishes. 

$1 = Ro(9)R 1 (6)<t>T c , 

A e = — {R 1 (e)d e R (e) + Ro(9)[Ri 1 (d)d e Ri(e)}R 1 (e)} , (2.19) 
igr 

where Ri{6) is chosen so that Ri(6) G Gi, i?i(0) = I and i? (27r)i?i(27r) = U 2 x d, with 
U2 G G2. All the fields are now single valued, and since Ri{6) annihilates &o ac , the form of 
the original part of the string solution is unaffected. 
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The same arguments apply to the subsequent symmetry breakings, leading to the ansatz 

= n^r, 

3=0 

2 n ft-1 "1 

^ = — e n^w^w^wi n . ( 2 - 2 °) 

^ i=0 [i=0 J=i-1 J 

where the gauge transformations i = 0...n — 1, satisfy G Gj, i?«(0) = I and 

U l j=0 R j( 2 ^) = U i+1 x d for some G 

If we define G n+ \ to be /, and choose an R n which satisfies the above conditions, the 

theory's gauge fields can be defined as 

= ft RMMQ) f[ R A e ) • ( 2 - 21 ) 

3=0 j=n 

Such a definition is single valued for all £ L(G\), and so the resulting cosmic string is 
not an Alice string. 

Although the conditions after (2.20) restrict the choice of Ri, they do not generally 
determine it uniquely. Thus it is necessary to look at the string field equations in detail. 

In the above discussion we have only considered the form of the string solution away 
from the string core. In some string solutions different parts of the Higgs fields have different 
radial and angular dependence (see section 2.3). In these cases the angular dependence of 
each part must be such that it is single valued. Also, since the solution must regular, all 
the winding parts of the Higgs fields must be zero at the string's centre. 

We will now apply the above arguments to the 50(10) cosmic strings discussed in 
section 2.3. These have i? o (0) = e inTsd ■ For all choices of T s , e 2™T s $vac = $ va C; go we can 
put R\{9) = /, which is the minimal energy choice. However e 2mnTs ^Jff / &\q c for some 
T s , so a non-trivial choice of i?2(#) may be required. To determine the most favourable 
choice we have to look at the form of the possible solutions close to the cosmic string. 

2.5 The Electroweak Symmetry Breaking 

Topological strings only form in a symmetry breaking G — ► H if tt\{G/H) ^ I. This is 
not the case at the electroweak ($io) symmetry breaking, so such strings do not form there. 
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It is still possible for <E>io to wind, and for string-like solutions to appear [36]. However, 
since it is energetically favourable for the Higgs field to unwind, they are not completely 
stable (although they could take a long time to decay). 

The situation is different in the presence of a topological string, formed at a previous 
symmetry breaking. As we discussed in section 2.4 the string gauge field may force a Higgs 
field to wind and take a string-like solution. Unlike the electroweak strings considered by 
Vachaspati [36], such solutions would be stable. 

As we discussed in the previous section, we need to add an extra singular gauge trans- 
formation to the cosmic string ansatz in order to give a single valued electroweak Higgs 
field (<J>io). We will define it to be e iTow6 '. Thus at large r 

$ 10 = e mT ^^ 6 ( Vu H u + Vd H d ) , (2.22) 

Aq = —T s + -e mT ^T ew e- mT ^. (2.23) 
gr gr 

Before we consider the radial dependence of the solution, we need to express <3?io and 
T ew in terms of eigenstates of the string generators (as we did with $126 in section 2.3). 

Thus we set Vu,d H u,d = £j,fc<%?> with T stfk = jtfk' T ew<%? = k<j>"£, and T cw = Y.j T j, 
with [T s ,Tj] = jTj. A suitable ansatz is then 

= E elin3+k)e (f jk h u jk {r) + 4>%h%{rj) , (2.24) 

A e = ^r s + y^e'"% je - inT " e . (2.25) 
gr j gr 

If [T S ,T CW ] = the second gauge term reduces to 26o(r)T cw /(gr). If [T s ,T ew ] 7^ the 
r component of the gauge field equations (2.6) becomes non-trivial and we also need to 
introduce a non-zero A r field, made up of e mTs jTje~ mTs terms. At r = the functions 
bj(r) and h u -^(r) take values which give a regular solution, and at r = 00 they are all equal 
to 1. Since $10 is single valued, nj + k must be an integer whenever cf^jf / 0. Any A r field 
must be zero at r = 00. 

We can now construct an approximate solution and use this to get an estimate of the 
energy. Defining r s to be the radius of the string (outside of which | <I> 126 1 takes its usual 
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VEV), and r cw to be the radius of the region in which |$io| does not take its usual VEV, the 
solution can be approximated separately in three regions. If all the coupling constants are 
of order 1, r s will be of order l^illl -1 = Vg 1 an d r cw will be of order l^ljfl -1 = Ve with 
Ve = Vu + Vd- ^10 takes the role of the Weinberg-Salam Higgs field and so rj E ~ 10 2 GeV. 

We will set all fields to their asymptotic values for r > r cw . For r s < r < r cw we set 
a(r) and $126 to their asymptotic values, and assume all other fields to be small. To first 
order the field equations are 

(hf k d )» + L£l _ k - h ^ = , h ] - A = . (2.26) 

These are solved by hf d {r) ~ r ± l fe l and bj{r) ~ r 2 or bj(r) = constant. For r < r s we take 
all fields to be small, and the field equations reduce to 

(h#y + ^p- - ^L^V = , - i = . (2.27) 

These are solved by h™ k (r) ~ r' n - J+fc ' and hj(r) ~ r 2 (the other solutions are not regular at 
r = 0). Requiring continuity of the solution and its first derivative at r = r s , reveals that 
the h U j k d (r) ~ r~l fc l and bj(r) = constant solutions can be neglected outside the string core. 
Requiring continuity at r = r cw gives the trial solution 



h% d (r) = ft fc (r) = { 



(r ew ) 



f / r \ 2 

r s < r < r c 



I r > 

&j ( r ) = 6 ( r ) = <j Uewy ■» - ■ — _ (2 2g) 
! I 1 r > r ew 



We assume any A r field is negligible. An estimate of the energy of this string-like solution 
can now be found by substituting the trial solution into the Lagrangian. The contribution 
from the region r < r s is suppressed by powers of r s /r cw ~ 1CT 14 and can be neglected. All 
the contributions are zero for r > r cw , thus 



E = 27r^°° rdr I |a r $io| 2 + |D^io| 2 + ^(a r ^) 2 + y| 

•£ w rdr Y, [k 2 + £(i - bfhi^ \r k + 4\ 2 + 



2vr / rdryZ\ h F £ + '-(I - bfHij \4> u k + ^ + \ (^) ' Tr(T 2 J + V 

n (\ k \ \ k \ 2k2 k2 ^lJL^i j.di2 

( 2 + 2 2 + 2\k\ + I+2\k\ J m + 41 



+ ^ T <Tl) + 2vr Vrdr . (2.29) 

y 'ew Jr s 
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where = £j cf>% d . 

The generator T cw should annihilate $^26 > since any change to the GUT string will give 
a large increase in energy. It is therefore a combination of electroweak generators. Gluon 
and photon generators do not affect &iq c , so the most energetically favourable choice will 
be of the form 

T cw = zr[z] + i= + w*t[w-}) , (2.30) 

where r[z] is the generator of the Z boson field, etc. z (which is real) and w (complex) are 
parameters to be determined. We will now find the eigenstates (4\' d ) of T cw which make 
up H® and H®, and the corresponding eigenvalues (A^' d ). 
Defining z' = \/iz, their eigenvalues are 



and the eigenstates are 



z' 2 + \w\ 2 



w 



(2.31) 



z ± ^Jz' 2 + \w\ 2 



z ± \ z' 2 + \w\ 



K + wH, 



u J 2^z rl + \w\ l 
Vd 



2^/z' 2 + \w\ 2 



Thus 



(2.32) 
(2.33) 



outside the string core. 

Substituting these results into (2.29), and taking gr cw /2 = r/" 1 gives 

E « 2^ ^(l A ±H^| 2 + l A ill4| 2 )+^ T e- 

= IOttt? 2 ,^ 2 + \ w \ 2 ) . (2.34) 

Thus the minimal energy solution will have z and \w\ as small as possible, subject to 
e mTs (f)± d being single valued. If 4>± d is made up of eigenstates of T s with eigenvalues [i± d 
then this is true if \± d + [i± d are all integers. 

If \w\ = the situation is slightly different. In this case ^jf is composed of 2 eigenstates 
of T ew , and not 4. They are H® and H%, and have eigenvalues 



\ u =°-z' , A d = -V 
5 5 



(2.35) 
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If H® and H® are made up of eigenstates of T s with eigenvalues /it"' , then it is only necessary 
for X u,d + to all be integers if $10 is to be single valued. 

2.6 Electroweak Cosmic String Solutions 

Using the results of section 2.5 we can now find the form of the electroweak Higgs and 
gauge fields around the various 50(10) cosmic string solutions. 

2.6.1 The Abelian U(l) String 

With the abelian string, T S H%~ = \H%- and T S H% + = -|#°* + . Thus e inT ° e <S>\f = 
r] u H®e me / 5 + r] l iH^ L e~ me ^ is not generally single valued, so a non-zero T cw is needed. 

The most general suitable T cw is given by (2.30). We can see that the eigenstates of 
T ew which make up (given by (2.32)) are also eigenstates of T s . The minimal energy 
choice of T ew must minimise Tr(T e 2 w ), while satisfying the conditions after (2.34) (or after 
(2.35) if \w\ =0). For this string the conditions are that 

A£ - n/5 and A± + n/5 (if w ^ 0) , (2.36) 
X d - n/5 and A" + n/5 (if w = 0) , (2.37) 

are all integers. This occurs when w = and ^8/5z+n/5 is an integer, or when 3z/ (2\/T0)± 
\/5z' 2 /8 + \w\' 2 + n/5 are integers. The choice which minimises (2.34) is w = 0, z = 
— \/5/8{n/5}, where {x} is defined as x minus the nearest integer. Thus if n is a multiple 
of 5, T cw = 0. 

Since [T S ,T CW ] = the electroweak string solution can be written in the form 
$ 10 = ( Vu Hy me h u (r) + Vd H d e- ime h d (r)) , 

A e = ^!lr s + ^T ew A,=0 otherwise, (2.38) 

gr gr 

where m = n/5 — {n/5}. The solution has the boundary conditions h u ^(oo) = b(oo) = 1, 
6(0) = 0, and if m 7^ then h u ^(0) = too. If m = 0, which occurs when |n| < 3, $10 need 
not be zero a the string's centre. We can get an estimate of it value using a trial solution 
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like that in section 2.5. The solution has h u ^{r) ~ (r/r ew )^ n//5 ^ for r s < r < r cw and 
/i M (r) constant for r < r s . By continuity at r = r B , ?i M (0) « (r s /r cw )^ n / 5 ^ ~ i(T 14 IW5}| _ 
Thus electroweak symmetry is almost fully restored for abelian strings even when $10 does 
not wind. Although symmetry is restored, $10 will only wind like a string if the GUT string 
has a high winding number (|n| > 2). This is in agreement with earlier work by Alford and 
Wilczek [37]. 

Taking $ 45 to be constant, the $10 terms of the potential will reduce to 

V = ^(|$ 10 | 2 -^) 2 + ^|$ 10 C$ 10 -2 w | 2 

+ A X (>i26| 2 -t7g) (l^iol 2 -^) (2.39) 

for the above string solution. C is a conjugation matrix which maps H" d to H2 U - In general 
the full potential will contain various cross terms such as |3>io • ^45 | 2 - With a suitable choice 
of parameters these terms will ensure that <3>io breaks SU(2)i rather than SU(3) C [13]. The 
resulting field equations are 



K,d [{fKl-fe)-f(l-a)] 



2 



K,d + ^ - U5JV ^ ^ 2 5 ^ - K A = Ax?? 2 (f - i)K, d 



+ (« + VdK - + A' E ^, u (/i«^ - l)h d , u , (2.40) 



6' 4 



(« + «) • (2-41) 



5^ 

The electroweak string will also contribute terms to the GUT string field equations. 
There will be an additional A x (?? 2 /i 2 + Vd^d ~ Ve)/ on the right hand side of (2.14) and a 
-g 2 [(1 - a)/5 - (1 - &){§ }/n] (r/ 2 ,/^ + term on the right hand gide of (2.15). These 

extra terms are far smaller than those already present in (2.14,2.15), so the back reaction 
from the electroweak string will be negligible. 

The field equations (2.40,2.41) can be solved numerically. The electroweak fields are 
shown in figure 2.2. The parameters used were %/r7 G = 10~ 14 - 5 , r) u = 3%, g = A E = A' E = 
A G = 1 and A x = 0. As expected, the effect of the electroweak fields on the GUT string 
was negligible. 

In general the size of the region of electroweak symmetry restoration is far greater than 
the GUT symmetry restoration, since r ew /r s ~ 10 14 . If n is a multiple of 5 then T ew = (as 
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h d {r) 

2 3 ^4 

Figure 2.2: Plot of electroweak field profiles for abelian string. 

would be expected). In this case |<3?io| will take its usual constant value for r > r s , and while 
the electroweak symmetry is still restored, it is over a much smaller region. When m = 
the electroweak Higgs field is non-zero at the string's centre. Its value can be evaluated 
numerically, and figure 2.4 shows its variation with respect to r/ G /ry E . It can be seen that 
h u ,d(0) ~ Vg/Ve when r/ G 3> r/ E , which is less than the value suggested by the trial solution. 

2.6.2 X s and X' SU(2) Strings 

Although at the GUT scale all the SU{2) strings have the same properties, this is not true 
at low temperatures, since they affect $io differently. The generators corresponding to the 
Xgj and X'^ both annihilate the usual vacuum expectation value of $>io, and so the string 
gauge fields have no effect on this symmetry breaking, and <I>io can be constant everywhere. 
There is still the possibility of symmetry restoration from potential terms, which is due to 
variation of the $126 and $45 fields (see section 2.7.3), although this is less significant. The 
other generators do have a non-trivial effect. 
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2.6.3 W R SU(2) Strings 

Different combinations of Wj^ generators give strings like those in section 2.3, but they 
are all gauge equivalent under U(1)q. Without loss of generality we can take T s = 
(t[w+] + t[h£~]) /(2\/2) and n > 0. For such a string, 

T s (H° uT H+)=±±(H° uT H+), 

T s (H° d ±H-) = ±±(H° d ±H~). (2.42) 

Thus e 2 ™ Ts cf>y(f = (-l) n $^ c , and so e inT * e $]ff is not single valued for odd n. As with 
the U(l) string we need a non-zero T ew . Again the most energetically favourable choice 
will minimise Tr(T c 2 w ), while ensuring that $io is single valued. If we split the eigenstates 
of T ew (2.32) up into eigenstates of T s , then $io will be single valued if the sum of the T s 
and T ew eigenvalues for each eigenstate are integers (see section 2.5). 

From (2.42) we can see that the T s eigenvalues are ±1/2. Thus we require that 

X^ d - n/2 and X u / + n/2 (if w + 0) , (2.43) 
X u ' d - n/2 and X u,d + n/2 (if w = 0) , (2.44) 

are all integers. The T ew eigenvalues are X± d or A"' d and are given by (2.31,2.35). These 
conditions are satisfied by T ew = if n is even (as expected). If n is odd then \± d + 1/2 (or 
just X u,d + 1/2 if w = 0) must be integers. The choice which minimises Tt(T^ w ) is z = 0, 
\w\ = 1/2. 

Using the trial solution (2.28) does not tell us what phase of w will give the lowest 
energy. We can find this by considering the contribution to the energy from the region 
r < r s . If \w\ = 1/2, the electroweak Higgs field is given by 

*io= E 4>ih jk (r)e^ 6 , (2.45) 
j,k=±i 



where 



-{^ - 2jkwn d ) [H° u - 2kw*H~ - 2jkw*H° d - jH+] . (2.46) 
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Matching solutions at r = r s gives hjk{r) ~ v /r s/ r cw(?*/^s)' riJI+fc ^ 2 - Its contribution to the 
energy in the string core is then 



E « 2vr / S rdr|a r $i | 2 + 
JO 



7T 



7T 



'"ew 



\n+ II 



10 



+ In - II 



(2.47) 



This is minimised by w = 1/2. For a more general choice of T s , this implies the gauge 
fields are related by W^/W^T = /W£~ . Combining the above results, and noting that 
[T S ,T CW ] = 0, leads to the ansatz 

1 



10 



J(nT B +T ew ) 



(H° u - H° d )r, + h + (r) + (H° u + H%)r,-h-(r) 



= (e im ' e (f>+ + e- im ' 6 Vl) /»+(r) + (e im Vt + e' mt V+) (r) , 



2no(r) 26(r) 



A M = otherwise 



(2.48) 



gr gr 

where i]± = (r] u =F rid)/V~2, m = (n — l)/2 and m! = (n + l)/2. For the solution to have the 
correct asymptotic form and to be regular at r = the required boundary conditions are 
h±(oo) = b(oo) = 1, h + (0) = 6(0) = 0, and also h-(0) = if n ^ 1. Ignoring any potential 
terms which couple <3?io to any of the other Higgs fields, the field equations for odd n are 



/4 [{l-b)±n(l-a)]\ 
h ± + ~ 4^2 h ± 



r 



^(vlhl + r?_tf_ - 4) ± ^ffaW - « + 2 %% ) 



* 2 | [(1 _ b)+ n(1 _ o)] ^ + [(1 _ fe) _ n(1 _ o)] ^2 | 



fc± , (2.49) 



(2.50) 



The field equations (2.49,2.50) can be solved numerically. The electroweak fields are 
shown in figure 2.3. The parameters used were the same as those for figure 2.2, thus 
rj + = \/2 x 10~ 14 ' 5 , n_ = 2r] + , and all coupling constants were set to 1. 

If n = 1, h-(0) need not be zero. As with the U(l) string we can estimate its value 
by extending the trial solution (2.28) to r < r s . This gives h-(0) ~ \/r s /r cw ~ 10~ 7 



so 



<I>io ~ inside the string when n = 1. Thus electroweak symmetry is fully or almost fully 
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2 4 rj E r 6 8 10 

Figure 2.3: Plot of electroweak field profiles for W^-string. 

restored there for all n. In contrast, the symmetry breaking caused by the GUT Higgs 
fields is only partially restored at the centre of the string, because |$i26| is about 7] G /\^2 
there. Figure 2.4 shows the numerically determined variation of h-(0) with respect to 
r] G /r] E . As with the abelian string, /i_(0) ~ Vg/Ve when r/ G S> i] E , which is less than the 
value suggested by the trial solution. 
If rj u = (2.48) simplifies, and 

$ 10 = e^ nTs+T ^ e ^ c h-(r) , (2.51) 

so $io has just one winding number, instead of the usual two. In this respect it resembles the 
equivalent abelian case. The field equations for h- and b are then the same as (2.49,2.50), 
but with 77-1- = 0. 

If n is even an electroweak gauge field is not required. Putting T ew = into (2.48) gives 
a suitable solution. In this case m = m! = n/2. The resulting field equations are then 
(2.49) with 6 = 1. In this case the region of electroweak symmetry restoration is limited to 
the string core. 

2.6.4 Y' SU(2) Strings 

The fourth type of nonabelian string has T s proportional to an appropriate linear com- 
bination of the generators of the gauge fields. The different combinations are all 
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Figure 2.4: Variation of electroweak Higgs field at the centres of the abelian and W R strings. 

equivalent under SU(3) C x U(1)q, and so without loss of generality we can take T s = 
(t[y{+] + t[y{-]) /(2>/2), and n > 0. As with the W^-string e 27rinTs ^ c = (-l) n $^ c , so we 
will need a non-zero T ew . Using 

T s (H° uT H l d ) = ±^(H° uT Hj), 

T s (H d ±Hl) = ± l -{H Q d ±Hl) , 
T S H-=T S H+ = 0, (2.52) 

we can split (f)± d into eigenstates of T s . This can then be used to determine the winding 
number of each component, and hence conditions for 3>io to be single valued (as we did 
with the other strings). For a general T cw with the required conditions are that 

A M + n ; ^ A M_^ (2 53) 

are all integers. If n is odd, no choice of T cw satisfies them. If however we take T cw to be 
proportional to the generator of the Z boson, then 4> u < d = H® d and there are no H~ or H d 
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components to worry about. This reduces the number of constraints. Now only \ u < d + t- 
and \ u,d — 7£ need to be integers. These conditions are satisfied by \z\ = \/T0/8. 

Unlike all the other cosmic strings we have considered, [T S ,T CW ] / for this string 
generator. This introduces several complications. As with the Higgs fields we have to 
split T ew up into eigenvalues of T 2 . If we define T\ = [T s , [T s ,T ew ]] and To = T cw — T\, 
then [T B , [T s ,Tj]] = Tj, as required. The angular dependence of Ag leads to a non-trivial r 
component of the gauge field equations (2.6). In order for them to be satisfied a non-zero 
A r field proportional to [T S ,T CW ] is needed. 

Additionally, we need to alter the angular dependence of $10 . Although introducing T ew 
terms gives a single valued $10 away from the string, it does not ensure that the eigenstates 
of T s 2 which make it up are single valued. If T cw = <jy/Wr[z]/8, with a = ±1, a suitable 
ansatz is 



* 10 = e^^ 9 e^Y,^i h u±(r) + ^h d± (r)) 



± 

gr gr gr 

A r = -i^lL inTse [T s ,T cw ]e- inTs(> , A (1 = otherwise, (2.54) 
9 3 

where m = (n — l)/2, m' = (n + l)/2, and 

< = f (AS T ^i) , ^ = f (flS ± ^i) • (2-55) 

Tg W is chosen to give a single valued <I>io inside the string, without altering the solution 
(2.20) outside the string. Thus it must be a combination of gluon and photon generators. 
Tg W = (t[g s ]/V3 + t[g 3 ] - t[a]/V6) /4 is a suitable choice. 

As with the other strings we expect the back-reaction of the electroweak string on the 
GUT string to be small. For simplicity we will neglect it. We will also ignore all the 
cross terms in the potential, and the effects of $45 (which couples to A r , but not the other 
electroweak fields). The field equations for the electroweak Higgs field and the additional 
gauge fields are 



2kh' uT +[k' + ^)h UT 
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[n(l-a)±(l-|6 )] 2 + (|6i) r 



vt(K + +K_) + r)j(hj + + hi. ) - 2T& h u ± 



+ 



4 
4 



K± ± l 
r 2 



K± + 



JL 

16r 2 



3(1 - |6o)6i 



2kh' dT + [k' + -)h dT 



- t k 2 + 



Kl-a)T(l-|6o)] 2 + (|fci) 2 ^ 



2/1,2 



4 
4 

+ y [K+hd- + h d+ h u _ - 2} r]lh UT , 



hd± + 



3(1 - f 6p)fti 
16r 2 



b 'i ,2, 8n 
r 3 



2a'k - (1 - a) (V - ^ 
2/o/i&o - (/o + /1 )&i + rilK+K- + r]jh d+ h d _ 



56 n 



K 'j fyf 2./O.A6, - (/* + /' ( 7)/; 



+ 10 
15 



2 

+10n 



?? 2 (/> 2 + + /* 2 _) + v 2 M + + ^ - 

r)lh u+ h u - + r]jh d+ h d _ b x 
ril(hl + -h 2 u _)- V 2 d (h 2 d+ -h 2 d _)} (1-a 



r^(l-a) 2 + ^ 



A; - — (61 a' + &;(!- a )) 



V <j ( V'cfo + 4 [«+ + K-) + rid(K + + 



2 2 > 

+ \(K + h' u _ - h u -ti u+ ) - ^(h d+ h' d _ - h d _h' d+ ) 
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The boundary conditions of the above functions are 6,(0) = h u+ (0) = hd-(0) = k(0) = 
fc(oo) = 0, bj(oo) = h u> d±(oo) = 1, and if n / 1 then /i M _(0) = /i<f+(0) = too. The 
potential terms considered in deriving the above equations are not enough to give the 
required boundary conditions at r = oo. The effect of $45 's coupling to <E>io also needs to 
be considered, although for simplicity we have neglected it in the derivation of the field 
equations. 

We find that the above equations are independent of a, the sign of T ew . The energy of 
the string is also independent, so the most favourable choice of electroweak gauge field is 
degenerate. As with the Wjj-string, if n is even no extra gauge fields are needed. The field 
equations are then the same as above, but with k = 0, hi = 1, h u+ = h u - and hd+ = hd- 
everywhere. The resulting electroweak symmetry restoration will be restricted to the string 
core. 

2.6.5 Summary 

Electroweak symmetry is restored and electroweak string gauge fields are present around 
the abelian and the W R and Y' nonabelian GUT strings. This generally occurs in a region 
around the string whose size is inversely proportional to the electroweak Higgs VEV, and 
is much bigger than the string core. If n is a multiple of 5 for the abelian, or 2 for the 
nonabelian strings, the region is approximately the same as the string core, and there are 
no extra string gauge fields. 

It is also possible that $10 will wind. For the abelian string its winding number is the 
closest integer to re/5, and hence zero for n = 1. For the W R and Y' nonabelian strings it 
is n/2 for even n. For odd n, different parts of $10 have different winding numbers, a bit 
like the corresponding GUT string. They are (n — l)/2 and (n + l)/2. The remaining two 
nonabelian strings (X s and X') have no effect on $10 at all. 

With the y'-string the electroweak gauge field has angular dependence, which leads to 
a non-zero radial gauge field component. Also, the choice of electroweak gauge field varies 
with distance from the string's centre. Both these additional effects are restricted to the 
string core. 
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2.7 Other Symmetry Restorations 
2.7.1 The Intermediate Symmetry Restoration 

So far the effect of the string on the second Higgs field $45 has been neglected, because it 
is far less significant. $45 is in the adjoint representation, and so its covariant derivative 
takes the form £^$45 = ^$45 — ^ig [A^, $45]. The generator corresponding to the Z' 
particle (P) commutes with $45, so the gauge fields of the abelian string will not stop $45 
from taking its usual vacuum expectation value everywhere. Thus it has no effect on the 
electroweak symmetry breaking, and gives no additional contribution to the energy. The 
other strings will give non-vanishing covariant derivatives at infinity. This is avoided by 
allowing $45 to wind like a string 

$ 45 = e ineTs ^{r)e~ indTs , (2.61) 

where $45 (00) is equal to the usual vacuum expectation value of $45. Conveniently e 2wniTs 
(for all choices of T s ) and commute, so <3?45 will be single valued for all n, and no extra 
gauge terms are needed. As with the <I>i26 and <E>io fields it is necessary to split $45 up into 
eigenstates of T s 2 . Thus, using the fact that $45° and T s 2 commute, and that T s 3 = jT s 

d> 45 = e ineT °ip lSl (r)e- ineT ° + ^ s (r) (2.62) 

with 

^ = 2T^\f - 2T S $™ C T S , (2.63) 

^0 = $ir " V>i • (2.64) 

The field equations for so,i( r ) will be similar to (2.11) and (2.12), with the similar boundary 
conditions so,i(°°) = 1 an d s\(0) = 0. Since the gauge contribution (1 — a(r)) vanishes 
for r > r B , $2 } (r) = $I| C is a solution outside the string. Thus the region of symmetry 
restoration will be of radius r s (order | ^* Y2 e I 1 ) - This is in contrast to the other Higgs fields, 
which restore symmetry in regions of order the reciprocal of their own values at r = 00. 

After the electroweak symmetry breaking the form of the solution for the Y'-string will 
be altered slightly due to the effect of the non-zero A r field. As with $126, we expect the 
effect to be tiny. 
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2.7.2 The Minimal Energy Choice of T s 

After the first symmetry breaking, there were only two gauge inequivalent possible strings, 
and the nonabelian strings had the lowest energy. After the second symmetry breaking there 
are 3 inequivalent types of nonabelian string. The Wjj-string has the highest energy, the 
X 5 -string has less, while the X' and Y' strings, which are gauge equivalent under SU(2)l, 
have the lowest energy. The final symmetry breaking gives an additional contribution to 
the y'-string, so the most energetically favourable choice of string generator is made up 
of rfx^], and does not have any effect on electroweak symmetry. Just because the other 
strings are not energetically favourable does not mean that they will not form, just that 
they are less likely to form (but see section 2.8). 

2.7.3 Non-gauge field symmetry restoration 

Even when a Higgs field is unaffected by a string's gauge fields, it is still possible for 
symmetry restoration to occur via the potential terms. This has previously been discussed 
for an abelian string in [30]. For example, before electroweak symmetry restoration occurs 
and in the absence of strings (so $io = 0), the potential takes the form 



where r] M is the usual VEV of $45, and |3>i26 ■ 3>45| 2 -like cross terms have been ignored. 
This is minimised by setting $126 an d ^45 to their usual VEVs. However, in the presence 
of an abelian string | ^> 126 1 is proportional to f(r), so writing |<3?45| = s(r)rj M , the potential 
becomes 



which for small r is no longer minimised by s = 1, if A' x / 0. Thus if the theory's parameters 
take appropriate values, | <E»45 1 will be lower than usual, or even zero at the string's centre. 
Alternatively |<E>45| could be higher than usual there. 



V($126, $45) = -J- (|$126 



2 \ 2 , A M / 1 |2 2 \ 2 
V G ) +^j-(J $ 4 5 | -Vm) 




(2.65) 




(2.66) 
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If symmetry restoration by this mechanism occurs at all, it will only be in the region 
r < r s , since |3>i26| and hence V(&u6, ^45) take their usual values at larger r. Unlike the 
corresponding symmetry restoration by gauge fields, $45 will not wind. If there is a very 
strong | < I > i26 • ^45 1 2 term present, <3?45 may be forced to wind. Such terms cause $45 to be 
orthogonal to $126) so it may be more favourable for $45 • $126 to remain zero than for 
$45 to have no angular dependence. When $10 is not equal to its usual VEV, it could also 
cause I $45 1 to vary, although this effect will be very small for most parameter ranges. In 
this case the symmetry restoration could take place in the larger r < r ew region. 

A similar situation can occur with <I>io in the presence of an X' or X s string. In this case 
both I $45 1 and |<l?i26| are lower than usual for r < r s . For nonabelian strings the potential 
is more complicated since it involves /1 and /o terms, as well as the corresponding $45 
terms. Even inside the string, /o and sq are non-zero, so the variation of the potential is 
likely to be less substantial than the abelian case, and hence extra symmetry restoration is 
less likely to occur. Even if it does, <J? 10 will take its usual VEV outside the string. Again 
we do not expect <I?io to wind, unless there is a very large |<&io • <&45 1 2 term, or something 
similar. 

2.8 Other Related Grand Unified Theories 

Although only one particular 50(10) GUT has been discussed, many of the results apply 
to different symmetry breakings. Any theory of the form 

50(10) ■ ■ ■ ^ . . . SU(3) C x SU(2) L x U(1) Y x Z 2 ^ SU(3) C x U(1) Q x Z 2 (2.67) 

could have string solutions of the form (2.7) or (2.8), which would cause electroweak sym- 
metry restoration at low temperatures in the same way as (2.1). The form of the other 
Higgs fields will not make much difference, as long as they are single valued in the presence 
of a string (like $45). If they are, it will not be necessary to add extra gauge fields, and so 
3>io will have the same behaviour as in (2.1). 

The Higgs fields which gain their VEVs after <E>i26 will determine the most energetically 
favourable choice of string generator, as $45 did in (2.1). If a GUT of the form (2.67) has 
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Higgs fields which take non-zero VEVs before $126, the choice of T s will be more restricted. 
If a generator has already been broken, the formation of the corresponding string will not 
occur. 

One theory of the form (2.67) is 

50(10) ^ SU(5) x U(1) P ^ SU(S) C x SU(2) L x U(1) Y x U(1) P 
^ SU{3) C x SU(2) L x U(1) Y x Z 2 

^ SU(3) C x U(1) Q x Z 2 . (2.68) 

The $a Higgs field transforms under either the 45 or 210 representation of 50(10), and is 
an SU(5) singlet. Unlike (2.1), only abelian strings can form in this theory, since the only 
generator that $126 breaks is P. This means that electroweak symmetry restoration will 
always occur in the presence of a string. <E>^ and $45 will both take their usual VEVs, so the 
only symmetries restored in the string core will be U(l)p, and the electroweak symmetry. 
Another interesting feature of this theory is that strings can form at energies close to the 
electroweak scale, and so r s could be of similar size to r ew , although still smaller. This string 
is a candidate for defect mediated electroweak baryogenesis [31]. Switching the second and 
third symmetry breakings also gives a theory with similar solutions. 

A different unifying gauge group (instead of 50(10)), with similar properties to (2.68) 
is SU(5) x U(l)p. It was suggested in ref. [11], and has two independent gauge coupling 
constants. Unlike (2.68), strings of all winding numbers will be topologically stable, since 
U(l)p is broken to Z instead of Z 2 . They could still decay by splitting into several strings 
with lower winding numbers. The field equations for the electroweak fields will be the same 
as (2.40) and (2.41), but with 1/5 the ratio of the two couplings instead of just 1/5. If the 
ratio is a, then <l?io's winding number will be the nearest integer to an/5 (with half integers 
rounded towards zero), so if |a| > |, $10 will always wind in the presence of a string. 

A theory which is substantially different from (2.1) starts with the symmetry breaking 
5O(10) -H SU(4) x SU{2) 2 x Z$ . The Z 2 C symmetry is not the Z 2 symmetry in (2.67). $i 
is not invariant under it, so it must be broken during or before the electroweak symmetry 
breaking. This will lead to formation of domain walls, and so such a theory will have 



43 



substantially different properties to (2.1), and is ruled out cosmologically [16]. 

Another type of theory closely related to (2.67) occurs when <I>i26 is replaced by <&i6, 
where the usual VEV of $i6 is proportional to ei. The gauge fields all gain masses in 
the same way as the equivalent theory involving $126 > but there will be no discrete Z2 
symmetry, so there will be no topological strings. However, solutions of the form (2.7,2.8) 
can still form, although since e 2mnTs will need to map ei to ei to give a single valued $16, 
only solutions with even n will occur. Of course, if such strings are to be observed, they 
will need to be stable, which will only happen for certain values of the theory's parameters. 
Embedded defects similar to these have been discussed previously [36, 38]. 

Yet another set of related theories can be obtained from (2.1) by choosing a different 
VEV of $45. Adding a multiple of P to it will not affect any of the SU(5) symmetry breaking 
since all the SU(5) fields commute with it, thus it will not alter which gauge bosons become 
superheavy. It will alter the sizes of the masses of the 50(10) fields. The most energetically 
favoured choice of nonabelian string will be the one with the lowest energy contribution at 
the $45 symmetry breaking. This will be the one whose string generator corresponds to 
gauge fields with the lowest mass. So by choosing $45 appropriately, a different nonabelian 
string could become most favourable. The X' and Y' strings are gauge equivalent at this 
stage, but since the X' string contributes nothing at the electroweak symmetry breaking, it 
will always be more favourable than the Y' string. Thus any of the X' , X s or W^-strings 
could be energetically favourable. If it is the Wfj-string, then it is most probable that 
electroweak symmetry restoration will occur. The same sort of freedom does not exist with 
<I>126 and $10, since any such change will give different fermion mass terms, and radically 
alter the theory. 

2.9 Conclusions 

In this chapter we have uncovered a very rich microstructure for 50(10) cosmic strings. In 
particular, we have found four nonabelian strings as well as one abelian string. We have 
examined the effect of the strings on the subsequent symmetry breakings. Our results arc 
summarised in the table 
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Gauge Type 
field 



Symmetry restoration 
50(10) 5*7(5) EW 



Z' 


U(l) 


yes 


no 


yes 




SU{2) 


partial 


partial 


no 




SU{2) 


partial 


partial 


no 




SU{2) 


partial 


partial 


yes 




SU(2) 


partial 


partial 


yes 



EW 
fields 



± 



W L 

Z (+ others) 



EW 
windings 



±(n/5 - {n/5}) 



±(n/2 ± {n/2}) 
±(n/2 ± {n/2}) 



We have defined {x} to be the fractional part of x, with |{x}| < 1/2. It seems that 
electroweak symmetry restoration by GUT strings is quite likely. The exact results are 
dependent on the details of the theory and the choice of string generator. For the 50(10) 
theory considered, electroweak symmetry is restored for the abelian string, and half the 
possible nonabelian strings, although the most energetically favourable of these does not 
restore electroweak symmetry. However, other closely related 50(10) GUTs have different 
minimal energy string solutions, which will restore the symmetry, such as (2.68), or (2.1) 
with a different choice of $45. Thus our results generalise to a range of theories. 

The size of the region of electroweak symmetry restoration for the topologically stable 
(n = 1) strings is determined by the electroweak scale, and is much larger than the string 
core. For nonabelian strings, with higher winding number, the region will be the same 
if they are topologically equivalent to the n = 1 string (i.e. odd n), and restricted to 
the string core if they are topologically equivalent to the vacuum (i.e. even n). There is 
no such distinction between topological and non-top ological abelian strings, which restore 
symmetry in the larger region if the winding number is not a multiple of 5. Some of the 
SU{5) symmetry is also restored by all of the nonabelian strings, but not the abelian string. 
This is only within the string core, irrespective of the string winding number, and since $45 
is not zero there, the restoration is only partial. 

For any choice of T s , ignoring possible potential driven symmetry restoration, the GUT 
will not be fully restored at the string's centre. All the nonabelian strings have non-zero 
(although smaller than usual) $126 and $45 fields at their centre, so the 50(10) (apart from 
the electroweak fields) symmetry is only partially restored inside the string. For the abelian 
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string <I>i26 is zero in the string core, but $45 takes its usual value, so with the exception 
of U(l)p, most of the 50(10) symmetry is broken. The resulting gauge boson masses are 
smaller, but still superheavy. However, for the abelian and the W R and Y' nonabelian 
strings, there is almost full restoration of electroweak symmetry in a larger region than the 
string core. 

It is also possible for additional symmetry to be restored by the potential terms of the 
theory. However this is far less significant than that arising from the string gauge field, and 
will require some tuning of parameters. 

Although the profile of the electroweak Higgs field obeys the same boundary conditions 
as a string, its exact form has a closer resemblance to a string with non-integer winding 
number. For the abelian string the actual winding number of $10 is less than that of the 
GUT string (about 1/5). The same is true for the nonabelian string, which has the winding 
number (or numbers if n is odd) of $10 about 1/2 that of the string itself. The winding of 
the electroweak Higgs field has implications for the existence of massless fermion currents, 
as we will discuss in chapter 3. 

In our analysis we have only considered terms occuring in the tree level Lagrangian. 
One-loop corrections are likely to induce couplings between the nonabelian string field and 
the electroweak Higgs. This may result in electroweak symmetry restoration around the 
X' and X s strings. However, the electroweak Higgs field is unlikely to wind in this region. 
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Chapter 3 

A Fermion Zero Mode Index 
Theorem 

3.1 Introduction 

It has been realised in the past few years that many cosmologically significant effects take 
place inside the core of strings [10]. One example is the formation of currents. These may 
provide a method of detecting strings. If they form at high energy scales, and are charged, 
the resulting electromagnetic field may be detectable [11]. It is also possible that currents 
will allow vortons (stable string loops) to exist. This may conflict with observations, in 
which case the corresponding theory can be ruled out [28]. There are several processes 
which can create currents on strings. These include interaction with the plasma, and 
collisions between cosmic strings. Charged currents can also be generated by magnetic or 
electric fields. 

In this chapter we investigate massless fermion currents. We do this by looking for 
non-trivial zero energy fermion solutions, or zero modes [26]. If they exist it is trivial to 
show that the string has light-like fermion currents (see section 1.6). 

As has been shown in chapter 2, subsequent phase transitions can have a considerable 
effect on the microphysics of cosmic strings. Zero modes on cosmic strings can be both 
created [32] and destroyed, thus creating or destroying the currents on the strings. If the 
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current is destroyed by subsequent micro-physical processes, then loops of cosmic string will 
no longer be stable, and vorton bounds on the theory will be evaded. In this chapter we 
discuss the fate of vortons and string superconductivity as the strings encounter subsequent 
phase transitions in a systematic fashion. 

In section 3.2 we derive an index theorem giving the number of zero modes for a general 
mass matrix. Whilst index theorems have been derived before [39, 40] they have been 
more restrictive in their validity and have only been able to determine the difference in 
right-moving and left-moving zero modes. Our index theorem has much more general 
applicability and can give a bound on the number of zero modes. 

In section 3.3 we discuss the existence of zero modes on all the different strings formed 
at the breaking of the 50(10) Grand Unified symmetry discussed in chapter 2. We show 
that at the electroweak phase transition zero modes can acquire a small mass which leads to 
dissipation of the string current. This allows vortons to decay and weakens the cosmological 
bounds on such models [41]. 

In section 3.4 we consider the implications of spectral flow. An important feature that 
allows zero modes to be removed is the presence of a massless particle (possibly another 
zero mode) that mixes with the zero mode after the transition. The implications of such 
couplings for current build up before the transition are also considered. 

We generalise some of the results of section 3.3 to other related theories in section 3.5. In 
other models the zero modes survive subsequent transitions allowing the associated vortons 
to persist. Such behaviour is displayed by a toy model discussed in section 3.6, where 
we explicitly construct the zero mode solutions after the symmetry breaking. Finally, we 
summarise our conclusions in section 3.7. 

3.2 Zero Mode Index 

Cosmic strings form in models with vacuum manifolds which are not simply connected. For 
example in a U(l) model, with potential (\(p\ 2 — i] 2 ) 2 , stable solutions exist with <f> = r\e %e 
at r = oo. In order for the total energy to be finite, a non-zero gauge field is needed to 
give a vanishing covariant derivative at r = oo. In a more general theory, involving a larger 
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group, G, string solutions take the form 

Hr,0) = e^U(r) , Ag = i-T(r) , (3.1) 

where T s is a generator of G that is broken by 4>. The choice of T s is restricted by the fact 
that (f> must be single valued. <j>(r) is equal to the usual VEV of 4> at r = oo, and must be 
regular at r = 0. T(r) obeys T(0) = 0, T(oo) = T s . 

In a general theory, T s will affect different components of (f> differently. This means that 
the various parts of <fi can have a wide range of winding numbers. In a theory with multiple 
phase transitions, the additional Higgs fields will be affected in the same way. It may also 
be necessary to alter T s at phase transitions to make the new Higgs fields single valued. 

In a theory with rif two-component fermions, the fermionic part of the Lagrangian is 

Acrmions = ^ a i^D^ a - -i^ a M a ^ + (ll. C.) , (3.2) 

where ip^ = ia 2 ^. If M a p depends on 9, as would be expected if M a p arose from the Higgs 
field of the string, then it is possible that the field equations will have non-trivial zero energy 
solutions. Solutions with only r and 9 dependence can be split up into eigenstates of a 3 : 
ipa • Such solutions have zero energy. If we solve the equations of motion in the background 
of a cosmic string, the field equations become 

e *° (d r + ify + eA( ^ + Ma/3l pL* = o ; (3.3) 

e- w (d r - l -d e - eAe^j ^ - M a ^f = , (3.4) 

where Ag is the string gauge field. If z and t dependence is added to the solutions they 
will correspond to currents flowing along the string. Their direction is left for those cor- 
responding to (3.3), and right for (3.4). In order to be physically relevant the solutions 
must be normalisable. Let Nl and Nr be the number of such solutions to (3.3) and (3.4) 
respectively. We attempt to derive an expression for them by generalising the analysis in 
ref. [26], which involves removing the 9 dependence of the problem, and then considering 
solutions near r = oo and r = 0. 

Choose the if) a s to be eigenstates of the string gauge field, with eigenvalues q a . The q a 
will depend on the fermion charges and the winding numbers of the various components of 
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the Higgs fields. Since the mass terms in (3.2) are gauge invariant, the angular dependence 
of the mass matrix must be 

M a p(r,6) = C a p{r)e iiqa+q P )e (no summation) . (3.5) 

The 9 dependence can also be factored out of the ip a s. 

^ = S^ e {U^{r)e m + Vt{r)e~ ue ) , (3.6) 
= j(i<*+\) e (U*(ry ie + V**{r)e- il9 ) . (3.7) 

I can take any value which gives a single valued ip. Thus I is integer or half-integer when 
the q a are half-integer or integer respectively. 

First consider left moving zero modes. Putting (3.5,3.7) into (3.3) gives equations for 
U% and V^. As r -> oo, C a/3 = 0(1) and eA e — ► T s /r, so 

i &r + Yr + I) V « + Ca ^ U P = • ( 3 - 9 ) 
Diagonalising C a p, we find that the 2nf solutions to these equations are modified Bessel 
functions. To leading order they are proportional to e ±Air /y / r, where Aj are C a ps eigenval- 
ues. Half of these large r solutions are normalisable at r = oo. If Aj = the corresponding 
solutions are r ±l ~ 1 / 2 . If I ^ (which is certain if the q a are not half integer), half of these 
solutions are acceptable. If I = the states corresponding to the solutions are physical, but 
are not localised to the string, and so are not of interest to us. 

Thus if all X { / or the q a are not half integer, exactly rif of the large r solutions 
are normalisable at r = oo. We will assume that this is the case in the following analysis. 
If these conditions do not hold, the index theorem derived below may overestimate the 
number of allowed solutions, and so will only provide an upper bound for the number of 
zero mode solutions. 

The Higgs fields, and hence M a p, are regular at the origin, so as r — > 0, C a/ 3 = 0(1), 
and eAg = 0(r). Thus 

[dr ~ qa ~} + l ^j + C a ^ L = , (3.10) 
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(d r - q -^-^j V a L + CaP Ui = . (3.11) 

To leading order, the small r solutions are 

JjL ^ r q a ~\+i f 

~ 0(l)r**+3+* V/3, (3-12) 

~ 0(l)r^+i+' V/3 / a , 

where each choice of a = 1 . . . rif gives one complex solution, and 

F £ _ r 9««-H ? 

U% ~ 0(l)r fo+ H V/3, (3.13) 
^ L ~ 0(l)r to+ H V/3 / a . 

This gives a total of 2ri{ independent complex solutions. For given I and a, (3.12) will be 
normalisable (for small r) if I > —q a + 1/2. If I < q a — 1/2 then (3.13) will be normalisable. 
Thus for a given I the number of well behaved small r solutions is 

™ f r i r n 

iV°(0 = ^ I I < q a - i + I Z > - 9a + \\ (3.14) 

a=l 

where I [X] equals 1 if X is true, and if X is false. What we are actually interested in is 
the number of solutions that are normalisable for all r (Nl(1)). Each such solution will be 
equal to some combination of the rtf well behaved solutions to (3.8,3.9) at large r, and a 
combination of the N®(1) suitable solutions to (3.10,3.11) for small r. If there are only nf, 
or less, suitable small r solutions, then in general any combination of the large r solutions 
will not be well behaved at r = 0. If there are nf + k suitable small r solutions, then k 
independent combinations of the large r solutions will be well behaved everywhere. It may 
be possible to get more solutions by fine tuning the theory, in which case the index derived 
would be a lower bound. 

The number of normalisable solutions for a given I is 

N L (l)=[N° L (l)-n i ] + , (3.15) 

where [x] + is defined to be equal to zero if x < 0, and x if x > 0. 
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This is not true if the equations obtained from (3.3) can be split into several independent 
sets. This will occur when M Q/ g is a direct sum of mass matrices. In this case the mass 
matrix can be split up into smaller matrices, which can be analysed individually. Even when 
M a n is not a direct sum of other matrices, it may still be possible to split the equations 
into two independent sets. This case will be considered separately later. 

Since U£ and are determined by real equations, each complex solution gives two real 
solutions. This suggests that the total number of left moving zero modes, Nl, is 2 J2i 
However, as can be seen from (3.7), solutions for I = k and I = —k are equal. For / = 
U% = ±V£, so one of ^'s solutions is zero. Thus the total number of independent real 
solutions is 



i i 



£,(l[l<q a -l]+l[l>-q a + i\)- nt 



a=l 



(3.16) 



+ 



The summation is over all values of I that give single valued ip. Since all the Higgs fields 
which make up M a p are single valued, (3.5) implies all q a or all q a — 1/2 are integers 
(assuming M a p is not a direct sum of smaller matrices), in which case respectively 1 — 1/2 
or I is an integer. 

A similar analysis can be applied to right moving zero modes. For large r the behaviour 
is the same. For small r, solutions are well behaved \i I > q a + 1/2 ot I < —q a — 1/2. This 
gives 



7lf 



< 



La=l 



+ 1 



I > 



rif 



(3.17) 



If q a is positive then one or two of the I [. . .] terms in (3.16) will be non-zero. If q a is 
negative, one or zero of them will be non-zero. By splitting q a into positive and negative 
eigenvalues, and ordering them, (3.16) and (3.17) can be simplified. If there are n+ positive 
and n_ negative q a s then, after reordering, q a = (g+, q£ . . . q£ + , q^ , q^ ... q^_, ... 0), 
where q^ > 0, q~ < 0. Clearly if the string gauge eigenvalues are not integer, there are no 
zeros. The I [. . .] terms in (3.16) and (3.17) can be combined to give 



Eih;+3<j<9;-3]-E^-3<i<-«7+3 



(3.18) 
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E 1 ^ + k<i<-oj -kl-^ht -\<i<^ + k 



(3.19) 



Li=i j'=i 
These expressions can be further simplified if the eigenvalues are ordered. If q± > q£ > 
■ ■ ■ > q_n + > an( i ?r — I2 — • • • — In- < 0, it is possible to evaluate the / summation by 
considering cancellation of the and qj terms. This gives 



min(n_ ,n+) 

*l = E 2 

min(n_ ,n+) 



E 2 

j'=i 



-«7 ~Qi 



+ E 2.;, 

j=n_+l 

+ - E 2,7 

j=n++l 



Taking the difference of these results gives 



(3.20) 
(3.21) 



J A = N L -N R 

min(n_,n+) n+ n_ 

- e %;+«7)+ e 2 ^;+ e 2 ?7 

J=l jr'=n_ + l j=n + + l 

n { -. 

= E 2 ^ = ^-[ lndetM ]^ • 

a=l 



(3.22) 



This is in agreement with other index theorems obtained elsewhere [39, 40]. The other 
index theorems were obtained by a different method and only gave 1a, not Nl and Nr. 
The total number of zero modes is 



min(n_,n+) 

1 = N L + N R = E 2 



n+ 



+ E 2 <T E 2^7, 

j=n_+l j=n++l 



(3.23) 



where only one of the last 2 terms contributes, depending on whether n + or n_ is bigger. 
This is also true of (3.22) and (3.20,3.21). 

If 1 is to be zero, then for every positive q a there must be one negative qp with the 
same magnitude. If every fermion field couples to a Higgs field with winding number zero, 
this will be the case. 

The above approach fails if C a p is of the form 

(3.24) 
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If C a p is assumed to have no zero eigenvalues, A a p and B a p are both nxn matrices, where 
n = nf/2. In this case when (3.5,3.7) are substituted into (3.3), two independent sets of 
equations are obtained. Expressions for Nl and Nr can found by considering just one set 
of these solutions. Putting 

U^(r)e m a = l...n 
a ' (3.25) 

V£*(r)e- M a = n+l...nt 

and (3.5) into (3.3) gives (3.8,3.9) for large r and (3.10,3.11) for small r, but with a more 
restricted range on the indices. In this case the allowed values of I need no longer be 
integers or half-integers, but the difference between any 2 values is still an integer. For 
large r, n of the nf complex solutions are normalisable. For a given I at small r, there is 
one normalisable solution for each q a satisfying I > —q a + 1/2 (a = 1 . . . n) or I < q a — 1/2 
(a = n + 1 . . . rif). 

Matching the solutions for large and small r gives 
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5Z(i['< 



Qa+n 



la=l 



+ 1 



I > -q a + \]) 



n 



(3.26) 



The solutions for different I are independent (unlike the previously considered cases), so 
the total number of real solutions (Nl) is just twice the number of complex solutions, thus 
Nl = 2J2i A similar expression can be obtained for Nr. 

The I [. . .] terms in (3.26) can be combined. Defining = —qj and qf = qj+ n for 
j = 1 . . . n, (3.26) and the corresponding expression for Nr, become 



n 

2(l[-^+i<Z<^-i]-l[^-i<Z<-^+i]) 



n 

2(i[^ + i<z<-^-i]-i[-^-Hi<^ + i]) 



+ 



(3.27) 



(3.28) 



If the string gauge eigenvalues are reordered, so that q^ > . . . > q£ and qf < . . . < q^, the 
expressions reduce to 



N L = £2[ g / + 



if 



(3.29) 
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N R = £2[-gf-gf] , (3-30) 

3=1 

which are similar to (3.20) and (3.21), with n + = n_ = n. This is not identical to the 
previous result, since the q a are divided up differently. 

When there are just two fermion fields involved, all the results reduce to 

J A = 2(gi + q 2 ) , (3.31) 
1 = |2" A | • (3.32) 

3.3 Fermion Zero Modes on SO(10) Cosmic Strings 

One example of a phenomenologically credible grand unified theory (GUT) has the sym- 
metry breaking 

50(10) ^ 5*7(5) x Z 2 (3.33) 

-H SU(3) C x SU(2) L x I7(l) y x Z 2 (3.34) 

^ SZ7(3) C x U(1) Q x Z 2 . (3.35) 

This theory was discussed in greater detail in chapter 2. Appendix A contains details 
of the fermion fields and their mass terms. A significant feature of its fermions are the 
massive neutrinos. It has been conjectured that rather than being massless, the left-handed 
neutrino actually has a very small mass. It would then provide a dark matter candidate. 
Recent experimental results suggest that neutrinos are massive, giving further credibility 
to the above model and those like it [42]. We can use these results to obtain an estimate 
of the right-handed neutrino mass and the 50(10) breaking scale. We will consider the 
top quark/tauon family, since renormalisation effects (which we will neglect) are likely to 
be less significant. The left-handed neutrino mass is approximately equal to m 2 jm v ( jR y 
The results of ref. [42] indicate that mixes with another type of neutrino which has 
a different mass, thus at least one of these neutrinos must be massive. Assuming the 
detected mass difference is the result of v T <-> mixing, Am 2 = m^r^ — tt?,^ = 
(5 x 10~ 4 — 6 x 10~ 3 )eV 2 . We will assume that m u ^ <C m„( r ), and so m 2 ,^^ ~ mf/ Am 2 . 
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Taking m t = (180 ± 12)GeV [43], we obtain m v(jK) « (4 - 16) x 10 14 GeV, suggesting that 
the ratio of the grand unification and electroweak scales is approximately (2 — 9) x 10 12 . 

We will now apply the results of the previous section to the various cosmic strings of 
this model. In several cases the string Higgs field (<3?i26) causes the electroweak Higgs field 
(3>io) to take a string-like solution. It is then possible that <l>io, like $126 w ih give rise to 
fermion zero modes [11, 32, 44]. The fermionic part of this theory's Lagrangian is 

Acrmions = ^L^D^ L - ^ E *L$10*£ - £*<7g* L$126* £ + (h. C.) , (3.36) 

where = (—J, a 1 ) since $i is a two component spinor. For simplicity we will only 
consider one family of fermions. Varying ^ l in (3.36) gives the field equations 

ia^D^L - ig E $ 10 ia 2 y* L - igo^iaHl = . (3.37) 

In section 3.2, an expression for the number of zero mode solutions to this type of equation 
was derived. The numbers and types of fermion zero modes can easily be found once (3.36) 
is put into the same form as (3.2). This involves calculating the mass matrix, M, for each 
type of string, and splitting it up into irreducible parts. Then the fermion fields coupling 
to each part, ip, need to be expressed in terms of eigenstates of T s . The eigenvalues of 
these states are then divided into two sets and inserted into the appropriate expressions 
(3.20,3.21,3.23) and (3.29,3.30). The way the eigenvalues are divided up depends on the 
mass matrix. Usually they will be split into positive and negative eigenvalues (see remarks 
before (3.20)), but if the mass matrix has the same type of degeneracy as (3.24), they are 
split according to which part of the matrix they couple to (see remarks before and after 
(3.27,3.28)). 

3.3.1 Zero Modes for the Abelian String 
High Temperature Neutrino Zero Modes 

At high temperatures $io is zero, and so with the exception of u°, none of the fermion fields 
are affected by Higgs fields. Thus the relevant part of the theory is just a two-component 
spinor coupling to an abelian string. This case has been discussed in section 1.6, and in 
greater detail in ref. [26]. 
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There are \n\ normalisable zero energy solutions. If z and t dependence are added, they 
all move to the left along the string if n > 0. If n < they move to the right. Thus conjugate 
neutrino zero modes always exist at high temperatures in the presence of an abelian string. 
For r > r s (assuming, for simplicity, that all coupling constants are approximately 1) the 
solutions decrease exponentially, so the zero modes are confined to the string core. 

High Temperature Non-Neutrino Zero Modes 

Although there is no Higgs field acting on the other fermion fields, it is possible for zero 
modes to be generated by the string gauge fields, as discussed by Stern and Yajnik [44]. 
The index theorem discussed in section 3.2 does not apply in this case, since it assumes 
that all the fermions considered couple to Higgs fields. 

Labelling the upper and lower components of the fermion fields X L and A^ respectively 
(as in (3.3,3.4)), where A = Ui,df, etc. (not u c ), the fermion field equations reduce to 



where a L,R = ±1, according to which component is being used. p\ is the eigenvalue of the 
field with respect to P/10: p e ~ = p v = pd\ = —3/10, p u <r = p Ui = Pdi = P e + = 1/10. There 
are normalisable solutions if \np\\ > 1, all of which can be found analytically. The number 
of solutions is equal to the highest integer that is less than Thus |n| must be at least 

4 for any zero modes of this type to exist. If the only stable strings have winding number 
1, then only conjugate neutrino zero modes will be present at high temperatures around an 
abelian string. 

Low Temperature Non-Neutrino Zero Modes 

At lower temperatures <I>io is non-zero and couples to all the fermion fields. Now none of 
the particles are massless, and all zero modes can be found using the index theorem. The 
electroweak phase transition changes the string generator used in section 3.2 to T s + T cw , 
with T ew defined in section 2.6.1. Before applying the theorem, the mass matrix needs to 
be split into irreducible parts. In this case there are 8 of them, one for each particle type. 




(3.38) 
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With the exception of the neutrino fields, there is just one Higgs field coupling to them. 
The Lagrangians for the non-neutrino fields are then (3.2), with 



where A can be di, Ui, or e~ . We have defined m u = <7 e t/ u /4, m e = ma = <7e%/4, and 
h e {f) = hd(r). The upper sign is taken for the itj and u\ fields, which couple to the H® 
component of $io- The lower sign applies for the di, d\ and e ± fields, since they couple to 
the H® component. Applying (3.32) reveals that there are 2\m\ solutions per particle type 
(14|m| total, not counting neutrinos). The field equations have previously been discussed 
in ref . [44] . 

If n is not a multiple of 5, the solutions decay exponentially outside r = r ew . When 
n is a multiple of 5 (in which case m = n/5), they decay outside r = r s . Thus the zero 
modes are confined to the region of symmetry restoration. The difference in sign between 
up and down quarks in (3.39) has physical significance when z and t dependence are added 
to the solutions. The up quark currents flow in the opposite direction to the down quark 
and electron currents. 

Low Temperature Neutrino Zero Modes 

The situation is more complex for the neutrino fields since they are affected by two Higgs 
fields at the same time. In this case (3.39) is replaced by 



where m G = g G rj G . This case is more complex than (3.39), although (3.32) still applies. It 
implies there are 2\m\ zero modes. Surprisingly this does not depend directly on n. As with 
the other particle zero modes, they will be confined to the region of symmetry restoration. 
They move left if m > (which only happens if n > 3), or right if m < (n < —3). 

Considering all fermions, there are a total of 16|m| zero modes, half of which are left 
moving, half are right moving. For a topologically stable string m = 0, so there are no 



M = 




(3.39) 




(3.40) 
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zero electroweak modes. In the case of the neutrinos this is slightly surprising, since at 
higher temperatures when $10 is zero, the abelian string does have neutrino zero modes. 
Intuitively, since $10 ~ {Ve/Vg)&126 ~ 10~ 14 $i26) the situation could be expected to be 
the same for lower temperatures. 

Since \2m\ < \n\, some of the neutrino zero modes will be destroyed by the electroweak 
phase transition. For a stable n = 1 string all zero modes are destroyed. Thus, since 
higher n strings usually decay, there are zero modes before, but not after the electroweak 
phase transition. It is expected that the neutral current in the string will disperse [45] 
in which case any vortons formed will dissipate after about lCT 10 sec [41]. Before the 
electroweak phase transition from about 10 10 GeV - 10 2 GeV the universe would undergo a 
period of matter domination. Once the vortons dissipate there would be some reheating of 
the universe. However the electroweak interactions and physics below the phase transition 
would be unaffected. 

3.3.2 Zero Modes of the SU(2) Strings 

There are two additional complications with nonabelian strings. Firstly the particle states 
are not eigenstates of the string generator, although this is easily solved by re-expressing the 
problem in terms of gauge eigenstates. Secondly, there are effectively twice as many Higgs 
fields, since each Higgs field has parts with two different winding numbers and different 
profiles. 

High Temperature Neutrino Zero Modes 

At high temperatures the gauge fields are proportional to T s . Since v c is not an eigenstate 
of T s , the problem must be expressed in terms of x c<y± ^ = iy c ± 2T s i/ c )/v / 2, which arc 



eigenvectors of T s . Their eigenvalues are ±1/2. The relevant part of the Lagrangian is 
(3.2), with 




(3.41) 
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Since /i(oo) = /o(oo) = 1 one of the fermion fields is massless at large r. This means that 
(3.32) does not apply in this case, as its derivation assumed that either the mass matrix has 
no zero eigenvalues at large r, or that T s does not have half integer eigenvalues. However, 
(3.32) can still be used to give an upper bound on the number of zero modes. Applying it 
to (3.41) reveals the upper bound to be zero. Hence there are no neutrino zero modes on 
the SU (2) strings for any values of n. 

High Temperature Non-Neutrino Zero Modes 

For the fermion fields that do not couple to $126 it is possible for zero modes to exist by 
the same mechanism as (3.38). However, unlike the abelian case, some fermion fields are 
annihilated by T s , so p\ is effectively zero, and they cannot have zero energy solutions for 
any value of n. For instance, the Ui, di, v and fields are all zero eigenstates of the string 
generator for the high temperature Wjj-string. Thus, in the presence of this type of string, 
solutions can only occur for the conjugate fields. Defining = (A ± 2T B X)/y/2, where 
T S A is not proportional to u c or T s u c , or equal to zero, the nonabelian equivalent of (3.38) 
is 



The solutions are normalisable if < I < ±na L ' R /2 — 1. Thus the total number of solutions, 
per number of particles (6 in this case), will be the largest integer below n/2. In order for 
any such solutions to exist n must be at least 3, so they do not occur for topologically 
stable strings. 

Low Temperature X s and X' String Zero Modes 

At low temperatures <E>io is non-zero and couples to all the fermion fields. When T s is made 
up of generators of the X s or X' fields, $10 just takes its usual vacuum expectation value. 
For the fermion fields that are not affected by $126 there is effectively no string and so no 




(3.42) 



It has the solutions 




(3.43) 
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zero modes. For the fields affected by $126 the solutions of the field equations will be at 
least as divergent as those of (3.41), so there will no normalisable solutions. 



Low Temperature W R -String Non-Neutrino Zero Modes 

The neutrino and electron fields all couple to $126 in the presence of a Wjj-string, while 
the quark fields are only affected by <&io. The string generator is now of T s + T cw . Its 
non-neutrino fermion eigenstates are x| ± ' > = (ui ± di)/y/2 and xt^ = ( u t T d%)/V%, where 
« = 1 ... 3. They are eigenstates of both T s and T ew , with T cw defined in section 2.6.3. The 
fermion mass matrix is obtained from (2.48). It is reducible into 4 parts. One couples the 
neutrino and electron fields, and the others couple the x| ± ' ) to the corresponding xt^ • The 
appropriate expressions to insert into (3.2) are then 

/ n n m , h , o im ' e m h P im ® \ 



M 
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(3.44) 



where m± = (m u =f m^jl. When n is odd the winding numbers of the Higgs components 
are m = (n — l)/2 and m! = (n + l)/2 

In this case the second version of the index theorem is needed since (3.45) is in the same 
form as (3.24). If m u / m^, (3.29,3.30) reveal that there are 2m = n — 1 left moving and 
2m right moving zero modes. If m u = m^, the e± im ' d terms are not present in (3.45), since 
m + = 0. The mass matrix is then reducible. Applying (3.32) to the 2 parts shows there 
are 2m left and right moving modes in this case as well. 

Since there are 3 choices of i there are a total of 12m different zero modes for the 
Wft-string after electroweak symmetry breaking (assuming no neutrino zero modes). The 
solutions are contained in the r < r e „ region. Since m = for the energetically stable n = 1 
string, it has no fermion zero modes. 

When n is even (so the string is actually topologically equivalent to the vacuum), M 
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takes the same form as above, but with m = m! = n/2. Thus the results for the odd n 
strings can be applied to even n strings, and there are 12(n/2) normalisable solutions. 



Low Temperature Y' String Non-Neutrino Zero Modes 

The y'-string can be approached in a similar way to the W^-string, although there are 
additional complications due to the form of the gauge fields. To put the problem in the 
same form as (3.2), it needs to be expressed in terms of the fermion eigenstates of the string 
generator, which in this case is T s + T cw + T' cw . T cw and T e ' w are defined in section 2.6.4. 
Defining Vi = (u\, u\, u\, v) and w\ = (e + , — c/3, di, —d1)a (where a is the sign of T ew + T^ W ), 



the eigenstates are Vi, Wi and \ 



(±) 



(v° ± wf) I a/2- Obtaining the fermion mass matrix 



from (2.54) we find, for i = 1 ... 3, that 
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, (3.45) 



This is similar to the mass matrix for the V^-string. However in this case the electroweak 
gauge fields couple to \\ ^ Although their effect is zero away from the string, they 
still mix the fermion solutions. Thus, unlike (3.45), the first form of the index theorem 
(3.23) is used. We find there are no normalisable zero mode solutions. If n is even there 
are no electroweak gauge fields, and the situation is the same as the W^-string. There are 
then 6n zero modes. 



Low Temperature W R and Y' String Neutrino Zero Modes 

For the fields affected by both $126 and $10 in the presence of a Wjj-string, the mass matrix 
is (3.45), with the M of (3.41) replacing the top left 2x2 submatrix. The eigenstates are 



now xl 

t} = {v± e~)/y/2 and X i = K =F e + )/\/2. This time however, the first version 



c(±) 
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of the index theorem (3.23) is used. It implies that there are no normalisable zero mode 
solutions for any choice of the parameters. 

For the Y'-string we take the fermion eigenstates to be V4, and X4~^ (defined above). 
The situation is similar to the W^-string, and so none of the SU(2) cosmic strings have 
low temperature zero modes involving the conjugate neutrino field. 

3.3.3 Summary 

The only fermion zero modes that form at high temperatures are v c zero modes around 
abelian strings (in which case there are |n| of them), or those that involve fermion fields 
that just couple to the string gauge fields, and not $126- This latter type of solution will 
only occur for higher n strings (n > 3 for SU(2) strings, |n| > 4 for abelian strings). 

At low temperatures there are a total of 16|m| different zero modes on an abelian string 
(\m\ for each particle type), where m is the winding number of $io- m = when \n\ < 3, so 
there are no zero modes around topologically stable abelian strings, and hence they can only 
be superconducting at low temperatures in the presence of an unusual Higgs potential [11]. 

If m / 0, and z and t dependence are added to the solutions, they will correspond 
to massless fermion currents. The electron and down quark currents will then flow in the 
opposite direction to the neutrino and up quark currents. 

In the presence of a X' or X s nonabelian string there are no zero modes at any tem- 
perature. After the electroweak phase transition the Wjj-string has 12m zero modes (m 
for each particle type not coupling to $126)) where m is the winding number of the part of 
§10 which winds least, m = n/2 for even n, and m = (n — l)/2 for odd n. The y'-string 
has same number of zero modes as the W^-string if n is even, and zero if n is odd. Thus 
for a minimal energy, topologically stable string there are no fermion zero modes, although 
there is still the possibility of superconductivity due to gauge boson zero modes. Thus even 
the X' and Xg— strings may be superconducting [46]. Indeed, it has been shown that such 
strings do become current carrying by gauge boson condensation [47]. 

The currents corresponding to any fermion zero modes present on the SU (2) strings 
do not consist of single particle types, as those around an abelian string do. Instead they 
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Figure 3.1: The Dirac spectrum of cosmic strings with a zero mode (left) and a very low 
lying bound state (right). Both spectra also have a bound state and continuum. 

are made up of eigenstates of the string generator. Also, unlike the abelian case, currents 
containing each particle type flow in both directions along the string. 

3.4 Index Theorems and Spectral Flow 

We have shown that zero modes can acquire masses at subsequent phase transitions. No 
matter how small this mass is, the spectrum of the Dirac operator changes significantly. 
If we compare the Dirac spectrum with a zero mode and a low lying bound state with 
infinitesimal mass (fig. 3.1), we see that an arbitrarily small perturbation to the zero mode 
introduces an entire new branch to the spectrum. Any massive state gives a spectrum 
that is symmetric about both the w and k axes, there is always a reference frame in which 
the particle is at rest and others where it is moving up or down the string. Conversely 
the zero mode, which is massless, can only move in one direction along the string and its 
spectrum is asymmetric. The transition from zero mode to low lying bound state causes 
drastic changes in the spectrum and can be brought about by infinitesimal changes in the 
value of one Higgs field. If we consider the species with the zero mode alone, this infinite 
susceptibility to the background fields appears unphysical. However, when we include the 
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massless neutrino in the 50(10) model the spectral changes are less worrying. For a small 
coupling between the two neutrinos, both the before and after spectra have a continuum 
of massless or nearly massless states. These states can be used to build the extra branch 
of the perturbed zero mode spectrum, allowing small changes in the overall spectrum for 
small changes in the background fields. 

This observation leads us to conjecture that zero modes can be removed only if they 
become mixed with other massless states. 

The coupling between the left and right handed neutrinos and the electroweak Higgs field 
need not be artificially small, the small mass of the light neutrino can be generated by the 
seesaw mechanism [48]. This coupling is present prior to the electroweak phase transition 
and allows transitions of the form vl + vr — ► //, where / is any light fermion from the 
standard model and the intermediate state is an electroweak Higgs. Such interactions allow 
zero modes on the string to scatter from massless neutrinos in the surrounding plasma and 
provide a current damping mechanism that affects current build up prior to the electroweak 
transition. 

3.5 Fermion Zero Modes in Other Theories 

As discussed in section 2.8, various symmetry breaking schemes give rise to the type of 
strings considered in section 3.3. The resulting strings will have the same kind of zero modes 
as theory (2.1), provided none of the other Higgs fields couple to the fermions. The only 
Higgs fields that can couple to fermions are those which transform under a representation 
contained in the 16 x 16, since fermion mass terms transform as a product of 16s. The 
only such representations are 126, 10, and 120 (which is antisymmetric), so the results of 
this section apply to a wide range of 50(10) theories. 

One alternative symmetry breaking is to use <E>i6 instead of $126 • Since $16 does not 
couple to the fermions, g G will be zero in (3.36), and so the neutrinos will have the same 
kind of zero modes as all the other particles (As would be the case if $126 were present, but 
g G were zero). However such a theory has left-handed neutrinos with significant masses, 
and an observable right-handed neutrino, so it is not compatible with the standard model 
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(unless some other mechanism is introduced to alter the neutrino masses). 

If the size of the coupling of the string gauge field to the electroweak Higgs field were 
different (as in the SU(5) x U(l)p theory mentioned in section 2.8), it is possible for \m\ to 
be non-zero on a topological string. In this case fermion zero modes can be present. This 
also means that neutrino zero modes may survive the electroweak phase transition. In fact 
if \2m/n\ > 1, extra neutrino zero modes will be created at the phase transition. 

Although we have concentrated on SO( 10), many of the results apply to a very wide 
range of theories. For example, we expect the destruction of zero modes at phase transitions 
will frequently occur. In a theory with Majorana mass terms, such as the right handed 
neutrino term in section 3.3.1, the first version of the index theorem can be used (3.23). 
After a phase transition, if every fermion couples to a non-winding Higgs field, such as $io, 
then all previously formed zero modes will be destroyed (see remarks after (3.23)). 

Another type of GUT with cosmic string solutions which couple to fermions are those 
with the symmetry group Eq [11, 49]. Such theories have a larger fermion sector (27 fields 
rather than 16), more phase transitions, and can have far more Yukawa couplings. Eq 
also has massive neutrinos, a total of 5 compared to .SC^lO^s 2. The details will be more 
complicated, but we expect currents to be formed and destroyed as in the 50(10) theories. 
We also expect the SU(2) strings of the theory to have zero, or very few, massless fermion 
currents. 

3.6 A Model with Persistent Zero Modes 

In this section we consider a model in which a zero mode survives a subsequent symmetry 
breaking, despite coupling to a Higgs field with a non-winding component. Consider a 
Majorana fermion similar to that in (1.16), but with a different Higgs field. In this case 
the Higgs field has two parts, a winding part from the string and a constant part from a 
second symmetry breaking, 




(3.46) 
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Such a Higgs field could only occur as the result of a global symmetry being broken. Note 
that because of this the above 4> cannot be written in the form used in (3.5), so the index 
theorem developed in section 3.2 does not apply in this case. We take only the upper 
component of tp to be non-zero, in which case its field equations reduce to 

(3.47) 
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(3.48) 
(3.49) 
(3.50) 
(3.51) 



We have explicitly constructed a zero mode after the second phase transition. The X 
coordinate is similar to the usual radial coordinate, but 'centres' on the effective zero of the 
resultant Higgs field, rather than the core of the string. 



3.7 Conclusions 

In this chapter we have seen that the microphysics of cosmic strings can be influenced by 
subsequent phase transitions. Fermion zero modes, and consequently superconductivity, of 
the strings can be created or destroyed by such phase transitions. In determining whether or 
not a cosmic string is superconducting it is not enough to just consider this at formation, 
but to follow the microphysics through the multiple phase transitions that the system 
undergoes. 

The existence of fermion zero modes at high temperatures enables the string to carry a 
current, and thus leads to the formation of vortons [24]. Normally, vortons formed at such 
high temperatures result in the theory being ruled out cosmologically [24, 28]. We have 
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shown that it is possible for currents formed at high energy to dissipate after a subsequent 
phase transition if the relevant fermion zero mode does not survive the phase transition. 
The vortons then cease to be stable, and cannot be used to rule out the theory [41]. This 
is the case with vr zero modes on an SO(W) string. Prior to dissipation there could be 
a period of vorton domination. After the phase transition the universe would reheat and 
then evolve as normal. 

To enable a systematic analysis of this effect we have derived a generalised index the- 
orem. Our index theorem is especially applicable to theories where the fermions acquire 
mass from more than one Higgs field. We applied the index theorem and also considered 
spectral flow. As a result we conjecture that zero modes are destroyed when they mix 
with other fermions that acquire mass at a subsequent phase transition from a non-winding 
Higgs field. 

We applied the index theorem to a realistic GUT, with an SO(10) symmetry group. 
The theory has 5 distinct types of string, all of which we considered. 

For the abelian string it is the winding number of <3?io that determines the existence of 
fermion zero modes after electroweak symmetry breaking. The number of zero modes is 
16 times its winding number, so there will be none for topologically stable strings, which 
have |n| = 1 and m = 0. Neutrino zero modes can always exist at high temperatures, but 
they do not survive the electroweak phase transition (for \n\ = 1). This has interesting 
implications as mentioned above. 

In the presence of a nonabelian string, different parts of <I>io can have different winding 
numbers. In the cases considered it is generally the part with the lowest winding number 
which determines the number of zero modes. Its winding number is equal to n/2 rounded 
down to the nearest whole integer for Wjj-strings, and for the Y'-strings with even n. The 
other strings have no zero modes. When present there are a total of 12 times this number 
of possible fermion zero modes. The fields coupling to $126 (part of which has winding 
number 0) do not have such solutions at any temperature. As with the abelian string there 
are no zero energy fermion solutions for topologically stable strings, and so fermion zero 
modes on strings are not as common as would be expected. 



68 



At high temperatures, it is also possible (for higher n) for fermion fields which couple to 
the string gauge fields (but not the Higgs fields) to have zero modes. The number of such 
zero modes is not determined by the index theorem as it only considers the effect of the 
Higgs field. However, such zero modes are easy to spot from the field equations. This can 
occur for the nonabelian strings, and for the non-conjugate neutrino fields around abelian 
strings. This effect is always overridden if the fermion field couples to a non-zero Higgs 
field. 

We expect qualitatively similar results in most other grand unified theories. 
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Chapter 4 

N — 1 Supersymmetric Abelian 
Cosmic Strings 

4.1 Introduction 

In this chapter, we investigate the microphysics of cosmic string solutions admitted by 
supersymmetric (SUSY) field theories [14]. This is important for at least two reasons. 
First, SUSY field theories include many popular candidate theories of physics above the 
electroweak scale. Second, the recent successes of duality in SUSY Yang-Mills theories 
may mean that the physics of nonperturbative solutions such as topological solitons may 
be easier to understand than in non-super symmetric theories. As in early studies of non- 
SUSY defects [8], we work in the context of the simplest models and in particular with 
versions of the abelian Higgs model obeying the supersymmetry algebra with one SUSY 
generator (iV = 1). We demonstrate that the particle content and interactions dictated by 
SUSY naturally give rise to cosmic string superconductivity in these models. Furthermore, 
by using SUSY transformations, we are able to find solutions for the fermion zero modes 
responsible for superconductivity in terms of the background string fields. A special case 
of the solutions discussed in this chapter has been obtained in a similar model by other 
authors using different techniques [50]. 

The effect of supersymmetry on the U(l)xU (1) bosonic superconducting Witten model 
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has been discussed by Morris [51]. The conclusion there is that simple bosonic supercon- 
ducting strings do not survive in the transition to the supersymmetric theory. Instead, in 
order to implement bosonic superconductivity in these theories, it is necessary to construct 
quite complicated models with a minimum of five chiral superfields. In the present chapter 
we study a different aspect of supersymmetric cosmic strings, namely the fermionic sector of 
the theories. In contrast to the results of ref. [51], we find that, even in the simplest SUSY 
abelian Higgs models, fermionic superconductivity is an inevitable result of the couplings 
and particle content required by the SUSY algebra. This powerful result leads us to the 
strong conclusion that all supersymmetric abelian cosmic strings are superconducting. 

The structure of the chapter is as follows. In section 4.2 we present the N = 1 SUSY 
abelian Higgs models. Such simple SUSY models are well-known in particle physics (for 
example see ref. [52]). However, we believe the cosmological relevance of the solutions 
we explore here to be new. In order to make contact with both the supersymmetry and 
cosmology literature, we employ both the superfield and component formalisms and repeat 
a number of well-established facts and conventions for the sake of clarity. Spontaneous 
symmetry breaking (SSB) in these models can be implemented in two distinct ways, leading 
to different theories with different particle content. We call these distinct models theory F 
and theory D respectively to refer to the origin of the SSB term in the Higgs potential. In 
section 4.3 we focus on theory F. We demonstrate how the cosmic string solution can be 
constructed in the bosonic sector, and determine the number of fermion zero modes using 
the index theorem in section 3.2. We then employ SUSY transformations to solve these 
equations in terms of the background string fields. In section 4.4 we repeat the analysis for 
theory D. The type of symmetry breaking in theory D is peculiar to theories with an abelian 
gauge group and we therefore expect theory F to be more representative of models with 
nonabelian gauge groups, such as grand unified theories. In section 4.5 we check our results 
for the special case discussed in ref. [50]. In fact, for theory D, the solutions are already of 
this special form. Finally, in section 4.6, we comment on the possible implications of our 
findings. 
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4.2 Supersymmetric Abelian Higgs Models and Spontaneous 
Symmetry Breaking 

Let us begin by denning our conventions. Throughout this chapter we use the Minkowski 
metric with signature —2, the antisymmetric 2-tensor £21 = e 12 = 1, 612 = e 21 = —1 and 
the Dirac gamma matrices in the representation 

\ 

(4-1) 





a» 



with (J M = (—I,a l ) and <r M = (—1,—a 1 ), where a 1 are the Pauli matrices. 

We consider supersymmetric versions of the spontaneously broken gauged U(l) abelian 
Higgs model. These models are related to or are simple extensions of those found in 
ref. [52]. In superfield notation, such a theory consists of a vector superfield V and m chiral 
superfields <E>j, (i = l...m), with U(l) charges In the Wess-Zumino gauge these may 
be expressed in component notation as 

V(x,0,0) = -{8^9)A^{x)+i9 2 8\{x)-i9 2 8\(x) + -9 2 9 2 D{x) , 
$>i{x,9,9) = fa(y) + V29^ i (y) + 9 2 F i (y) , (4.2) 

where y M = x M + i9a fl 9 [14]. Here, fa are complex scalar fields and is a vector field. 
These correspond to the familiar bosonic fields of the abelian Higgs model. The fermions 
tpia, X a and A a are Weyl spinors and the complex bosonic fields, Fi, and real bosonic field, 
D, are auxiliary fields. Finally, 9 and 8 are anticommuting superspace coordinates. In the 
component formulation of the theory one eliminates Fi and D via their equations of motion 
and performs a Grassmann integration over 9 and 9. Now define 

W a = ~D 2 D a V , (4.3) 
where D a and are the supersymmetric covariant derivatives and W a is the field strength 
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chiral superfield. The superspace Lagrangian density for the theory is then given by 

C = l (w a W a \ 2 + W & W% 2 ) + ^ie^ v ^\ d2§2 + W(*i)\e2 + W{*i)\g2 + ZD . (4.4) 

In this expression W is the superpotential, a holomorphic function of the chiral superfields 
(i.e. a function of $j only and not and W\g2 indicates the 6 2 component of W. The term 
linear in D is known as the Fayet-Iliopoulos term [53]. Such a term can only be present in 
a J7(l) theory, since it is not invariant under more general gauge transformations. 
For a renormalizable theory, the most general superpotential is 

1 1 

with the constants bij, Cij k symmetric in their indices. This can be written in component 
form as 

W(4>i,il>j,F k )\ e 2 = aiFi + bij (^Fi<f>j - 7;4>i4>j^ + c ijk {Fi4>j^ k - ipiipj4> k ) (4.6) 

and the Lagrangian (4.4) can then be expanded in Wess-Zumino gauge in terms of its 
component fields using (4.2). The equations of motion for the auxiliary fields are 

F* +a,i + bij<j)j + Cij k <j)j<j) k = , (4.7) 

D + £ + = . (4.8) 

Using these to eliminate Fj and D we obtain the Lagrangian density in component form as 

C = Cb + C f + C y -U , (4.9) 



with 



L B = (D^^D^-^F^, (4.10) 

C F = -i^D^i-iXa^d^X , (4.11) 

C Y = ^=q i (j)i'ipiX - (^-bij + Cijk^>kj ipiipj + (c. c.) , (4.12) 

U = \F { \ 2 + ^D 2 = \a i + b i j(t>j+c i jk<t>j<l)k\ 2 + , (4.13) 
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where D^ = + ^igqiA^ and = d^A u - d u A^. 

Now consider spontaneous symmetry breaking in these theories. Each term in the 
superpotential must be gauge invariant. This implies that ai 7^ only if qi = 0, bij 7^ 
only if qi + qj = 0, and 7^ only if qi + qj + q\~ = 0. The situation is a little 
more complicated than in non-SUSY theories, since anomaly cancellation in SUSY theories 
implies the existence of more than one chiral superfield (and hence Higgs field). In order 
to break the gauge symmetry, one may either induce SSB through an appropriate choice 
of superpotential, or, in the case of the U(l) gauge group, one may rely on a non-zero 
Fayet-Iliopoulos term. 

We shall refer to the theory with superpotential SSB (and, for simplicity, zero Fayet- 
Iliopoulos term) as theory F and the theory with SSB due to a non-zero Fayet-Iliopoulos 
term as theory D. Since the implementation of SSB in theory F can be repeated for more 
general gauge groups, we expect that this theory will be more representative of general 
defect-forming theories than theory D for which the mechanism of SSB is specific to the 
U{\) gauge group. 

4.3 Theory F: Vanishing Fayet-Iliopoulos Term 

The simplest model with vanishing Fayet-Iliopoulos term (£ = 0) and spontaneously broken 
gauge symmetry contains three chiral superfields. It is not possible to construct such a 
model with fewer superfields which does not either leave the gauge symmetry unbroken or 
possess a gauge anomaly. The fields are two charged fields <£±, with respective U(l) charges 
q± = ±1, and a neutral field, $o- A suitable superpotential is then 

W(*i) = c$ ($+$- - r] 2 ) , (4.14) 

with 77 and c real. The potential U is minimised when F{ = and D = 0. This occurs when 
4>o = 0, 4>+4>~ = T] 2 , and \4>+\ 2 = \4>~ | 2 . Thus we may write (f>± = r]e ±ta , where a is some 
function. It is interesting to note that U also has a local minimum at (f> + = </>_ = 0. <p$ 
is undetermined. This minimum has non-zero vacuum energy, and will cause the universe 
to expand exponentially. This process is called inflation [54], and is a potential solution to 
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many cosmological problems. It will stop as <f>o approaches 0. 

We shall now seek the Nielsen-Olesen [18] solution corresponding to an infinite straight 
cosmic string. We proceed in the same manner as for non-supersymmetric theories. Con- 
sider only the bosonic fields (i.e. set the fermions to zero) and in cylindrical polar coordinates 
(r, ip, z) write 

00 = 0, (4.15) 

0+ = </>*- = ve inip f(r) , (4.16) 

A, = --n^%, (4.17) 
g r p 

F± = D = , (4.18) 

F = C77 2 (l-/(r) 2 ) , (4.19) 

so that the z-axis is the axis of symmetry of the defect. The profile functions, f(r) and 
a(r), obey 

f+ f -- n 2 ^^ = cW( f - l)f , (4-20) 



a' 



a"-^ = -gW(l-a)f , (4.21) 

with boundary conditions /(0) = a(0) = and f(oo) = a(oo) = 1. Note here, in passing, 
an interesting aspect of topological defects in SUSY theories. The ground state of the 
theory is supersymmetric but spontaneously breaks the gauge symmetry, while in the core 
of the defect the gauge symmetry is restored but, since \Fi\ 2 / in the core, SUSY is 
spontaneously broken there. 

We have constructed a cosmic string solution in the bosonic sector of the theory. Now 
consider the fermionic sector. With the choice of superpotential (4.14) the component form 
of the Yukawa couplings becomes 

L Y = (M + - <^_V>-) A - c {M+il)- + <p+^- + <f>-^+) + (c c.) (4.22) 

As with a non-supersymmetric theory, non-trivial zero energy fermion solutions can exist 
around the string. The number of solutions can be determined with the index theorem 
derived in section 3.2. z and t dependence can easily be added to the solutions. They 
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then correspond to massless currents flowing along the string. The index theorem also 
determines their direction of travel. When the string solution (4.19) is substituted into 
(4.22), the non-zero terms can be written in the same form as (3.2), with 



M 













-ce ilp -ij-e* 













77/ (r) , V 



/ 



V A i 



(4.23) 



V *^ 

This mass matrix is in the same form as (3.24), so the second form of the theorem (3.29,3.30) 
is used. It shows that there are 2\n\ complex (or 4|n| real) solutions. Half their correspond- 
ing currents move left along the string, and the other half move right. 

In general, in non-super symmetric theories, it is difficult to find solutions for fermion 
zero modes in string backgrounds. However, in the supersymmetric case, SUSY transfor- 
mations relate the fermionic components of the superfields to the bosonic ones and we 
may use this to obtain the fermion solutions in terms of the background string fields. A 
SUSY transformation is implemented by the operator G = e € Q +€ Q, where e a are Grassmann 
parameters and Q a are the generators of the SUSY algebra which we may represent by 

d 



Qa — 



Q & = ^--ia^ad,*. 



(4.24) 
(4.25) 



In general such a transformation will induce a change of gauge. It is then necessary to 
perform an additional gauge transformation to return to the Wess-Zumino gauge in order to 
easily interpret the solutions. For an abelian theory, supersymmetric gauge transformations 
are of the form 



V 



e - iAqi <f>i , 

V + - (A — A) , 
9 



(4.26) 
(4.27) 
(4.28) 



where A is some chiral superfield. 
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Consider performing an infinitesimal SUSY transformation on (4.19), using d^A^ = 0. 
The appropriate A to return to Wess-Zumino gauge is 



A = igeo^OA^y) 



The component fields then transform in the following way 



(4.29) 



<f>±(y) 

2 F (y) 
-Oa^OA^x) 



cf) ± (y) + 2ie^eD^ ± (y) , 
e 2 F (y)+29eF (y) , 



(4.30) 
(4.31) 

-Oa^A^x) + ^ie 2 6a^a u eF^(x) - ^i&Oa^o 1 'eF^(x) . (4.32) 



Writing everything in terms of the background string fields, only the fermion fields are 
affected to first order by the transformation. They are given by 

2na' 



(V'o)a 

where we have defined 



gr 



^2 (if'a r T -(1 - a)faA e & V e ±inifi , 
V2c V 2 (l-f 2 )e a , 



(4.33) 

(4.34) 
(4.35) 



ie*v 



(4.36) 



-up 



e"? 



Let us choose e a so that only one component is non-zero. Taking ei 
—iS/(y/2r]), where 5 is a complex constant, the fermions become 



Ai = 5- 



grj r 



(V'o)i = 
ty_)i = 6* 



/' + -(!-«)/ 
r 



-iScn(i - f 2 ) , 



J(n-l)<p 



(4.37) 
and ei = 

(4.38) 

(4.39) 
(4.40) 
(4.41) 
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It is these fermion solutions which are responsible for the string superconductivity. Similar 
expressions can be found when e\ = 0. It is clear from these results that the string is not 
invariant under supersymmetry, and therefore breaks it. However, since f'(r),a'(r), l — a(r) 
and 1 — / 2 (r) are all approximately zero outside of the string core, the SUSY breaking and 
the zero modes are confined to the string. We note that this method gives us two zero mode 
solutions. Thus, for a winding number one string, we obtain the full spectrum, whereas for 
strings of higher winding number, only a partial spectrum is obtained. 

4.4 Theory D: Non-vanishing Fayet-Iliopoulos Term 

Now consider theory D in which there is just one primary charged chiral superfield involved 
in the symmetry breaking and a non-zero Fayet-Iliopoulos term. In order to avoid gauge 
anomalies, the model must contain other charged superfields. These are coupled to the 
primary superfield through terms in the superpotential such that the expectation values of 
the secondary chiral superfields are dynamically zero. One simple way to do this is add 8 
charge 1/2 chiral fields, X a , and a superpotential term: 



The extra fermions cancel the anomaly, and the scalar potential makes the bosonic parts 
of X a zero [55]. The secondary superfields have no effect on SSB and are invariant under 
SUSY transformations. Therefore for the rest of this section we shall concentrate on the 
primary chiral superfield which mediates the gauge symmetry breaking. 



Y,9a<5>X a X, 



a 



(4.42) 



a 



Defining £ 



r^gr] 2 , the theory is spontaneously broken and there exists a string 



solution obtained from the ansatz 



<P = Ve m *f(r), 



(4.43) 



A 




(4.44) 



(4.45) 



F = 



. 



(4.46) 
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Taking n > 0, the profile functions f(r) and a(r) obey the first order equations 

f = n { -^f , (4.47) 

4 = ^V(1-/ 2 )- (4-48) 

Now consider the fermionic sector of the theory, and apply the index theorem (3.32). For 
n > there are 2n zero modes, which move left. Performing a SUSY transformation, again 
using A as the gauge function to return to Wess-Zumino gauge, gives to first order 

A a — ±gv 2 (l-f 2 W + ° z tep, (4-49) 
^ — V2^(l - a)f(ia r - a^^e^ . (4.50) 

If ei = both these expressions are zero. The same is true of all higher order terms, 
and so the string is invariant under the corresponding transformation. For other e, taking 
ei = —iS/r] gives 

Ai = 6g V (l-f), (4.51) 
V>i = 2^*^(1 - a)/e i(n - 1)v . (4.52) 

Similar results are obtained when n < 0. In this case the zero modes move right, and the 
non-zero fields are A2 and ip2- 

In this theory supersymmetry is only half broken inside the string. This is in contrast 
to theory F which fully breaks supersymmetry in the string core. The theories also differ in 
that theory D's zero modes will only travel in one direction, while the zero modes of theory 
F (which has twice as many) travel in both directions. In both theories the zero modes and 
SUSY breaking are confined to the string core. 



4.5 The Super-Bogomolnyi Limit 

In non-supersymmetric theories it is usually difficult to find solutions for fermion zero modes 
on cosmic string backgrounds. In such theories one can, however, often obtain solutions in 
the Bogomolnyi limit which, in our theory, corresponds to choosing 
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2c 2 = g 2 . (4.53) 

In this limit, the energy of the vortex saturates a topological bound, there are no static 
forces between vortices and the equations of motion for the string fields reduce to a pair of 
coupled first order differential equations. It is a useful check of the solutions obtained in the 
previous sections to confirm that they reduce to those already known in the Bogomolnyi 
limit. 

Imposing (4.53) equations (4.20,4.21), together with the requirement of finite energy, 
become 

/' = n^(l-a), (4.54) 
a' 



n- = c Y(l-r) . (4.55) 



r 



Note that these are similar to (4.47,4.48). We see that all solutions to theory D are auto- 
matically Bogomolnyi solutions. Imposing (4.53) on (4.41) gives the following solutions. 



Ai = 5c V (l - f 2 ) , (4.56) 

(V>+)i = 25*nUl - a)^"" 1 ^ , (4.57) 
r 

(Vo)i = ~iSc V (l - f 2 ) , (4.58) 

(V>-)i = 0. (4.59) 

This limit, with n = 1, was considered for a similar theory by Garriga and Vachaspati [50] 
and the above results are in agreement with theirs. This is a useful check of the techniques 
we use. 



4.6 Concluding Remarks 

We have investigated the structure of cosmic string solutions to supersymmetric abelian 
Higgs models. For completeness we have analysed two models, differing by their method of 
spontaneous symmetry breaking. However, we expect theory F to be more representative 
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of general defect forming theories, since the SSB employed there is not specific to abelian 
gauge groups. 

We have shown that although SUSY remains unbroken outside the string, it is broken 
in the string core (in contrast to the gauge symmetry which is restored there). In theory F 
supersymmetry is broken completely in the string core by a non-zero .F-term, while in theory 
D supersymmetry is partially broken by a non-zero D-term. We have demonstrated that, 
due to the particle content and couplings dictated by SUSY, the cosmic string solutions 
to both theories are superconducting in the Witten sense. We believe this to be quite a 
powerful result, that all supersymmetric abelian cosmic strings are superconducting due to 
fermion zero modes. Note that this result is in contrast to those obtained in earlier analyses 
of purely bosonic superconductivity in Witten-type SUSY models [51]. An immediate and 
important application of the results of the present chapter is that SUSY GUTs which break 
to the standard model and yield abelian cosmic strings (such as some breaking schemes 
of 50(10)) may face strong constraints from cosmology [28]. The existence of zero modes 
around SUSY monopoles has been investigated previously [56]. However, in this chapter 
we have considered cosmic strings since, unlike monopoles, supercurrents on strings are 
cosmologically significant. 

Although explicitly solving for such zero modes in the case of non-supersymmetric 
theories is difficult, in the models we study it is possible to use SUSY transformations 
to relate the functional form of the fermionic solutions to those of the background string 
fields, which are well-studied. For theory D the solutions all obey the Bogomolnyi equations 
exactly, and for theory F we have also checked that the solutions we find reduce to those 
already known in the special case of the Bogomolnyi limit. 

While we have performed this first analysis for the toy model of an abelian string, we 
expect the techniques to be quite general and in fact to be more useful in nonabelian theories 
for which the equations for the fermion zero modes are significantly more complicated. The 
question of superconductivity in nonabelian SUSY cosmic strings is considered in the next 
chapter. We also examine the effect of soft SUSY breaking in all of the models. 
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Chapter 5 

Cosmic Strings and 
Supersymmetry Breaking 



5.1 Introduction 

In recent years, supersymmetry (SUSY) has become increasingly favoured as the theoret- 
ical structure underlying fundamental particle interactions. In this light it is natural to 
investigate possible cosmological implications of SUSY. 

In chapter 4 we discussed the effect of SUSY on the microphysics of simple cosmic string 
solutions of abelian field theories. In particular we developed and applied the technique 
of SUSY transformations to investigate the form of fermionic zero modes, required by 
SUSY, which lead to cosmic string superconductivity. In the present work we extend our 
original ideas to a more general class of field theories, namely those with a nonabelian 
gauge group. Since nonabelian gauge theories underlie modern particle physics and, in 
particular, unified field theories, this class of theories is a realistic toy model for grand 
unified theories (GUTs). The particular example we examine, SU{2) x U(l) — > U(l) x Z2, 
admits two types of string solution, one abelian and the other nonabelian. This model has a 
similar structure to 50(10) and should provide insight into cosmic strings in SUSY GUTs. 
Most of the features exhibited by this theory will also appear in larger nonabelian theories. 
We apply the technique of SUSY transformations to the nonabelian case and extract the 
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behaviour of the zero modes as functions of the background string fields. We then compare 
the results to those obtained in chapter 4 for abelian strings. 

Since SUSY is clearly broken in the universe today, it is important to know how the 
SUSY zero modes behave when soft-SUSY breaking occurs. We investigate this for both 
the abelian and nonabelian strings by explicitly introducing soft-SUSY breaking terms into 
the theory. The result is that all the zero modes are destroyed in almost all the theories, 
the exception being when a Fayet-Iliopoulos term is used to break the gauge symmetry in 
an abelian model. We briefly comment on the physical reasons for this and show how the 
effect may be seen through the breakdown of an appropriate index theorem. 

These results have a cosmological significance since fermion zero modes on the string can 
be excited, causing a current to flow along the string. The string then behaves as a perfect 
conductor. The existence of charge carriers changes the cosmology of cosmic strings. In 
particular, they can stabilise string loops, resulting in the production of vortons [24]. Such 
vortons can dominate the energy density of the Universe, and have been used to constrain 
GUT models with current-carrying strings [28]. However, if the zero modes are destroyed 
at the SUSY breaking energy scale, then the current, and hence vortons, can dissipate. 
Thus, the underlying theory may be cosmologically viable. 

The plan of this paper is as follows. In section 5.2 we construct a simple super symmetric 
model based on the group SU (2) x U(l) and display the abelian and nonabelian string 
solutions. In section 5.3, we use an index theorem to find its fermion zero modes. We 
then use SUSY transformations to obtain the zero modes in terms of the background string 
fields. Soft SUSY breaking terms are introduced in section 5.4, and their effect on the 
string solution is considered in section 5.5. In section 5.6 the effect of SUSY-breaking on 
the zero modes is analysed using the index theorem. In this section we consider both string 
solutions for the SU(2) x U(l) model and also for the U(l) theories discussed in chapter 4. 
For the nonabelian string the SUSY breaking terms destroy the zero modes, while for the 
other string solutions the situation is more complicated. 
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5.2 An SU(2) x U(l) Model 

There exist many nonabelian theories with breaking schemes giving rise to cosmic strings. 
In general both abelian and nonabelian strings can be produced in such a process, depending 
on which part of the vacuum manifold is involved in the winding. 

In this section we consider a simple example in which the gauge group SU(2) x U(l) is 
spontaneously broken down to the group £7(1) x Z 2 via the superpotential 

W = ci5 (* • * - rj 1 ) + c 2 {S$ T A$ + S$ T Aif>) . (5.1) 

The chiral superfields &i(4>i, ipi, Fi) and &i((fti,^i, Fi) are SU(2) triplets with U(l) charges 
±1 respectively. The other chiral superfields, S*o(so, too, F$ ), S(s,uj,Fs) and S(s, Cj, Fg), 
are SU(2) scalars with U(\) charges 0, —2 and +2 respectively. Finally, defining T 4 = 
^2/31, the vector supermultiplets are V a (A^, X a , D a ), a = 1,...,4. Since the constant 
matrix A satisfies AT* = -(T l )*A (i = 1,...,3), and the SU(2) gauge transformations 
are 5$ = iT a n a <& and <5<& = —iT a *n a &, the superpotential is gauge invariant. The scalar 
potential, derived in the standard manner [14], is then 

u = c 2 \cf> ■ 4> - v 2 \ 2 + - 4\ 2 + 4\24>i4>3 - 4> 2 2 \ 2 

e 2 

+ \c lS0 (f) + 2c 2 sA0| 2 + |cis o + 2c 2 sA0| 2 + — (|0| 2 - |0| 2 - 2|s| 2 + 2|s| 2 ) 2 

+ e 2 |(0! + cp 3 )<p* - (0x + 3 )&| 2 + y m 2 - \H 2 ~ \4>i? + \h\ 2 ) 2 ■ (5.2) 

This is minimised when all fields are zero except 4>i = (f>i = i] or at any (broken) gauge 
transformation of this. We note also that the theory has a local minimum with 4> = (f) = 
and that this structure can give rise to hybrid inflation [54] . This is true even for the abelian 
theory described in chapter 5. In both cases inflation ends with defect formation. 

As we mentioned above, there are abelian and nonabelian string solutions to this theory. 
The abelian solution is obtained from the ansatz 

>h = fa = Vf(r)e^ , (5.3) 

^ - • <«> 

F So = c lV 2 (l-f(r) 2 ) , (5.5) 
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where T G = sj\ T 3 + y'fT 4 . All other fields are zero and the profile functions a and / obey 
the boundary conditions /(0) = a(0) = and /(oo) = a(oo) = 1. 



The nonabelian solution is obtained from the ansatz 



A 




(5.6) 



(5.7) 



er 



F s 



Fs 




l Cl r, 2 (2-f(r) 2 -g(r) 2 ) , 



(5.8) 



(5.9) 



where e& are unit vectors obeying T 1 e/ C = ke^- In this case g(0) is finite, 5(00) = 1 and / 
and a obey the same boundary conditions as in the abelian case. 

Note that /, g and a are solutions to simple coupled second order ordinary differential 
equations. Their forms can be obtained numerically and are well known [34]. 

5.3 SUSY Transformations and Zero Modes 

The string solutions obtained above have all the fermion fields set to zero. In this section 
we investigate what happens when these fields are excited in the background of the cosmic 
string. We can find the number of zero modes with the index theorem in chapter 3. 

We already know that there must exist fermion zero modes on the string. Rather than 
attempting to solve the difficult fermion equations of motion to obtain them, we shall use 
the technique described in section 4.3 which exploits the power of SUSY to obtain some of 
the solutions. 

5.3.1 Abelian string 

The relevant part of the Lagrangian is the Yukawa sector which is entirely determined by 
supersymmetry. In the abelian case, the non-zero Yukawa couplings are 




+ 2c 2 (e^3W + e-^V3^) 
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r)f(r) + (c. c.) , 



(5.10) 



where X± = (A 1 =F iA 2 ) / a/2- With respect to the string generator, the only fields with 
non-zero eigenvalues are ipi and A+ (eigenvalue 1) and V>i an d A_ (eigenvalue —1). The 
Yukawa Lagrangian can be split up into 5 independent parts. Applying the index theorem 
(3.29,3.30) to these reveals that there are a total of three left moving and three right moving 
complex zero modes. 

Now, following the techniques in section 4.3, we perform an infinitesimal SUSY trans- 
formation (with Grassmann parameter e a ) followed by a gauge transformation to return to 
the Wess-Zumino gauge. The string fields all transform quadratically and so are unchanged 
to first order. However, the fermions transform linearly and, in terms of the background 
string fields, it is possible to find two complex (or 4 real) fermion zero mode solutions given 
by 



lj = V2c lV 2 (l- f 2 )e , 
X G = 



/3a 
-—a e , 
5 er 



V>i = iV2 V e^ ( f'a r + i- — -fa ( 



€ . 



(5.11) 
(5.12) 

(5.13) 

(5.14) 



Setting either component of e to zero gives one of the zero modes. One is left moving, the 
other moves right. 



5.3.2 Nonabelian string 

In the nonabelian case it is convenient to split tf>, ip and A* into eigenvectors of T 1 , and 
label them by their eigenvalues. Defining Xi = hPi + i>i)/V2, Q = (ipi — ^j)/\/2 and 
= (u> ± lj)/\/2, the non-zero Yukawa couplings are 



£y = -Cir/ 

- c 2 r] 

- c 2 r] 



^o-^(x + ^-X-e^)/(r)A c 



-V^Co5(r) + (C + e-^ + C-e^)/(r) 



W(_) 
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ert 

+ y 

ier) 



V2(C-A+ - ( + \-)g(r) + Co(A + e"^ - A"e^)/(r) 
V2( g(r) + (C+e-^ + C-^)/(r)l A 4 + (c. c.) . 



(5.15) 



In this case x±> C± an d A 1 * 1 have eigenvalues ±1. Applying (3.29,3.30) to the two irreducible 
parts of (5.15) shows that there are just two complex zero modes, moving in opposite 
directions. 

Once again performing an infinitesimal SUSY transformation and a (nonabelian) gauge 
transformation we obtain two complex zero modes 



^0 = -^c lV 2 (2-g 2 -f 2 )e, 

-.2/2 



W(+) = c 2 v {g - f )e , 
A = 
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-i — a e , 
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X+ = *Ve 



X- = ive~ ilp f'a r -i 



Xo = iV2rjg'a r e . 



r 

l-a 



/ ^ ) I , 



(5.16) 
(5.17) 
(5.18) 

(5.19) 

(5.20) 
(5.21) 



Thus in this case there are no zero modes beyond those implied by SUSY, in contrast to 
the abelian case. This is related to the fact that there are components of the Higgs fields 
that do not wind in the nonabelian case. 



5.4 Soft SUSY Breaking 

Perhaps the most attractive feature of supersymmetry arises from the non-renormalisation 
theorems, which provide a solution to the hierarchy problem. The mass of the electroweak 
Higgs field receives quadratically divergent radiative corrections from the particles which it 
couples to. In general this leads to conflict with experiment. In a supersymmetric theory 
the contributions from the fermion and boson fields cancel each other exactly. This ensures 
that quadratic divergences are absent, and so any tree-level hierarchy of scales is protected 
from receiving quantum corrections. 
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Super symmetry is not observed at everyday temperatures, and so must have been broken 
as the universe cooled. It is crucial that the quadratic divergences remain absent from the 
theory, so that the hierarchy problem is still avoided. This is achieved by adding only 'soft' 
SUSY breaking terms. These are all either mass terms, or couplings with positive mass 
dimension [15]. Obviously they must be non-invariant under SUSY too. All the allowed 
terms are equivalent to one of three types. They lead to the following changes to the 
Lagrangian 

1. Arbitrary mass terms for scalar particles are added to the scalar potential. 

2. Bilinear and trilinear scalar terms in the superpotential, plus their Hermitian conju- 
gates, are added to the scalar potential with arbitrary coupling. 

3. Mass terms for the gauginos are added to the Lagrangian density. 
For a general theory the allowed terms are 



-£ soft = myth + (tiijQiQj + dijkQ&jQk + (c c.)) + m x X a X a + (c. c.) . (5.22) 



Of course, only those soft SUSY breaking terms which respect the gauge and other sym- 
metries of the theory are allowed. Because of this the bilinear and trilinear terms resemble 
those of the superpotential, and are conventionally written as multiples of them. Experi- 
ment suggests that the SUSY breaking scale is around 1 TeV. 

5.5 SUSY breaking and an Abelian String Model 

In section 4.3 we referred to an abelian theory in which the gauge symmetry is broken via 
an F term as 'theory F\ The corresponding superpotential was 



W = c$ ($+$- - V 2 ) ■ 



(5.23) 



The allowed soft SUSY breaking contributions to the potential are 



mf \4>i\ 2 + cA4>o4>+4>- + (c. c.) . 



(5.24) 
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Defining m 2 = (m 2 + + m_)/2 and £ = {m + — m_)/g, the scalar potential becomes 

V = c 2 |^--r/ 2 | 2 + |0o| 2 (c 2 (|^-| 2 + |^ + | 2 )+m2)+m 2 (|^ + | 2 + |^_| 2 ) 

+ ^^ + | 2 -|^| 2 + ^) +c^o^- + c^(0o0 + ^-)*-y • (5.25) 

We see that £ acts like a Fayet-Iliopoulos term. Any stationary points of the potential will 
have 

CM ^- (5.26) 



m2 + C 2(|^ + |2 + |0_|2) • 

Thus (j)Q will acquire a non-zero expectation value if the trilinear term is present. 

It is not generally possible to find the minimum of (5.25) analytically. Instead we will 
consider a couple of special cases. It is convienient to define m = m/(cr)), £ = £,/{gr] 2 ) and 
P = 2c 2 /g 2 . 

If A = the potential has stationary points at (pi = and 



771 2 

6 ± = r?e^/l-- 

cosh x 

^0 = 0, (5.27) 
where a is real and % satisfies 

sinhx + m 2 (/?- l)tanhx + |= (5.28) 

and 

coshx-m 2 >0. (5.29) 

If m 2 > 1 and |£//3| 2 < (m 4 — 1) then (5.28) has no real solutions which also satisfy (5.29). 
The only minimum of V is then </>j = 0, and so the SUSY-breaking terms restore the broken 
U(l) gauge symmetry. 

If m 2 < 0, m 2 {[3 - 1) < -1 and ||| 2 / 3 < [-m 2 (/3 - l)] 2 / 3 - 1 then (5.28) has 3 solutions, 
two of which are minima. Since the potential has disconnected minima, it is possible for 
domain walls to form. Domain wall formation occurs when a discrete vacuum symmetry 
is broken. In this case it is the 4> + <-> <p*_ symmetry. If £ = the minima have the same 
energy. The domain walls are then stable, and will come to dominate the energy density 
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of the universe as it evolves, conflicting with observations [16]. If £ / the minima will 
have different energies, and the walls can decay by quantum tunnelling. If this happened 
rapidly enough, conflict with observations could be avoided. 

Other values of the parameters give a unique global minimum with <f>± ^ 0. Since we 
require the phase transition not to be reversed, and that no domain walls form, the range 
of allowed parameters is restricted. However, we expect m,( < 1 in general, so these 
restrictions are not very significant. If the phase transition and supersymmetry breaking 
occur close together then m,£ ~ 1, and the allowed parameter range may be significantly 
reduced. 

If £ = TOq = 0, the vacuum is given by 



If A = and ttiq < the potential (5.25) resembles that of the bosonic superconductivity 
model in ref. [11]. Although c/)q = outside of a cosmic string, it may be energetically 
favourable for it to be non-zero inside. Excitations of this condensate could then form 
currents flowing along the string, although these are likely to be small. 

5.6 Fermion Zero Modes after SUSY breaking 

We will now consider the effect of soft SUSY breaking on the fermion zero modes. Since 
the techniques we have used for finding the zero modes are strictly valid only when SUSY is 
exact, it is necessary to investigate the effect of these soft terms on the zero modes we have 
identified. As we have already commented, the existence of the zero modes can be seen 
as a consequence of the index theorem in section 3.2. The index is insensitive to the size 
and exact form of the Yukawa couplings, as long as they are regular for small r, and tend 
to a constant at large r. In fact, the existence of zero modes relies only on the existence 
of the appropriate Yukawa couplings and that they have the correct ip dependence. Thus 
there can only be a change in the number of zero modes if the soft breaking terms induce 




^ = 'Ye' 



(5.30) 



90 



specific new Yukawa couplings in the theory and it is this that we must check for. Further, 
it was conjectured in section 3.4 that the destruction of a zero mode occurs only when the 
relevant fermion mixes with another massless fermion. 

5.6.1 U(l) Abelian models 
Theory F 

We saw in section 5.5 that although the scalar mass terms will alter the values of <p + and 
(f>—, they do not generally produce any new Yukawa terms. Thus these soft SUSY-breaking 
terms have no effect on the existence of the zero modes. The possible exception was a 
negative nig term. This may lead to a <fo condensate inside the string which will produce 
a Y>+Y>- term. 

The presence of the trilinear term gives (fro a non-zero expectation value everywhere, 
and hence produces a Yukawa term coupling the Y>+ an d V'- fields. This destroys all the 
zero modes in the theory since the left and right moving zero modes mix. 

For completeness note that a gaugino mass term also mixes the left and right zero 
modes, aiding in their destruction. 

In terms of the index theorem, the change in the number of zero modes arises because 
(3.20,3.21) applies after the SUSY breaking, while (3.29,3.30) applied before it. Although 
the fermion eigenvalues do not change, the expression relating them to the zero modes does. 

Theory D 

The £7(1) theory with gauge symmetry broken via a Fayet-Iliopoulos term and no super- 
potential is simpler to analyse. New Higgs mass terms have no effect on the number of 
zero modes, as in the above case, and there are no trilinear terms. Further, although the 
gaugino mass terms also affect the form of the zero mode solutions, they do not affect their 
existence, and so in theory D the zero modes remain even after SUSY breaking. 
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5.6.2 SU(2) x U(l) model 
Abelian strings 

The effect of soft SUSY breaking terms on the zero modes which were found analytically in 
(5.14) is identical to the equivalent U(l) theory. Thus, all SUSY zero modes are destroyed in 
this case. In this larger theory there are also other, non-SUSY, zero modes. Not all of these 
are destroyed by a gaugino mass term, as some do not involve the gaugino fields. However, 
if the trilinear terms (or possibly a negative t«q term) give so a non-zero expectation value, 
all the zero modes are destroyed. The extra Yukawa terms mean there are fewer irreducible 
parts to the fermion mass matrices. This results in more terms cancelling in (3.29,3.30), 
reducing the number of zero modes. As in the other cases, the physical reason behind the 
destruction of the zero modes is that left and right movers mix. 

Nonabelian strings 

As in the other cases above, non-zero gaugino mass or trilinear terms destroy the zero 
modes that were found with SUSY transformations (5.21). Similarly (3.20,3.21) are required 
instead of (3.29,3.30), implying that the left and right moving modes mix. For nonabelian 
strings these are the only zero modes and so none remain after SUSY breaking. 

5.7 Comments and Conclusions 

We have examined the microphysics of abelian and nonabelian cosmic string solutions to 
the 50(10) inspired supersymmetric SU(2) x U{\) model. By performing infinitesimal 
SUSY transformations on the background string fields we have obtained the form of the 
fermionic zero modes responsible for cosmic string conductivity. These solutions may be 
compared to those found in chapter 4. 

Our results mean that fermion zero modes are always present around cosmic strings 
in SUSY. We conjecture that in theories with F-term gauge symmetry breaking, the zero 
modes given by SUSY always occur in pairs, one left and one right moving. It also seems 
likely that such theories always have hybrid inflation. 
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Furthermore, in the abelian case there were additional zero modes that were not a 
consequence of super symmetry. We expect that similar extra zero modes will be present in 
a larger theory, even in the nonabelian case. 

We have also analysed the effect of soft SUSY breaking on the existence of fermionic 
zero modes. The SU (2) x U{\) model and two simple abelian models were examined. In all 
cases Higgs mass terms did not affect the existence of the zero modes. In the theories with 
F-term symmetry breaking, gaugino mass terms destroyed all zero modes which involved 
gauginos and trilinear terms created extra Yukawa couplings which destroyed all the zero 
modes present. 

All the theories with F-term symmetry breaking feature a chargeless scalar field. If the 
SUSY-breaking contribution to its mass term is negative, a chargeless condensate may form 
inside the string. This would allow bosonic currents to flow along the string. Additionally 
the condensate would destroy the fermion zero modes, in a similar way to any trilinear 
terms. 

In the abelian theory with D-term symmetry breaking, the zero modes were unaffected 
by the SUSY breaking terms. It was conjectured in section 3.4 that zero modes would only 
disappear when they mixed with another massless fermion field and this is consistent with 
the results obtained in this chapter. If the remaining zero modes survive subsequent phase 
transitions, then stable vortons could result. Such vortons would dominate the energy 
density of the universe, rendering the underlying GUT cosmologically problematic. 

Therefore, although SUSY breaking may alleviate the cosmological disasters faced by 
superconducting cosmic strings [28], there are classes of string solution for which zero 
modes remain even after SUSY breaking. It remains to analyse all the phase transitions 
undergone by specific SUSY GUT models to see whether or not fermion zero modes survive 
down to the present time. If the zero modes do not survive SUSY breaking, the universe 
could experience a period of vorton domination beforehand, and then reheat and evolve as 
normal afterwards. 

If the zero modes do occur in pairs (one left and one right moving) in F-term gauge 
symmetry breaking, it is possible that they could scatter off each other [57]. This would 



93 



cause the current to decay, and could stop vorton domination. 

There is the possibility that even if zero modes are destroyed they become low-lying 
bound states. Such bound states may still be able to carry a persistent current. If this is the 
case, even such theories may not be safe cosmologically. Work on this is under investigation 
(and see chapter 6). 

It may also be possible to extend our analysis of the effect of SUSY breaking on the 
bosonic fields. The full potential is very complex, even in the abelian case. It may be 
worth using some sort of approximation or a numerical solution. We have already seen that 
domain walls can form in some special cases, allowing the corresponding parameter range 
to be ruled out. Other parameter ranges could be ruled out in the same way. 



94 



Chapter 6 



Massive Fermion Bound States 

6.1 Introduction 

In section 1.6 and chapter 3 we investigated fermion zero modes and massless currents 
on cosmic strings. It is reasonable to expect that cosmic strings may also have massive 
bound states and currents. It is possible that such currents occur in models without zero 
modes, in which case a greater range of models will have conducting strings. We expect 
the maximum size of these currents to be smaller than the massless case, since they will 
require less momentum to escape from the string. 

If the theory contains particles which are massless off the string, and these interact 
with the current carriers, we expect them to scatter into those states and thus the current 
will dissipate. This is most likely to happen in a GUT at high temperatures, since all the 
Standard Model fields are massless then. This is not the case after the electroweak phase 
transition, suggesting that massive currents will be most significant at low temperatures. 

In chapter 4 we showed that fermion zero modes occur in all supersymmetric cosmic 
string theories. It seems that they will gain a mass when supersymmetry is broken (see 
chapter 5), suggesting that massive fermion currents will occur in models with broken 
sup er symmetry. 

The existence of bound states will also affect scattering off the string. A stable bound 
state may help to catalyse certain interactions. This could provide a way to probe a string's 
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internal structure, and determine the properties of the GUT in which it formed. 

We found that each of massless fermion currents considered in chapter 3 only moved in 
one direction along a cosmic string. If a massive fermion current exists on a string, there 
will be equal numbers of left and right moving states. If lower energy states are always 
filled before higher energy ones, there will be equal amounts of left and right movers. The 
total current is then zero. The fermion states will still have non-zero angular momentum, 
and so they may still stabilise string loops to form vortons. A net current will develop if 
the states are not filled evenly, due to difficulties in dissipating momentum. 

In section 6.2 we will derive the field equations for fermion bound states in a general 
theory, and use them to show that such states are always time or light like. We solve these 
equations numerically for an abelian model in section 6.3. In section 6.4 we speculate about 
the existence and type of bound states in some other theories. The results are summarised 
in section 6.5. 



6.2 Field Equations and a Positivity Condition 

Consider a general fermion Lagrangian, 

Acrmions = ^ai^D^a ~ ^^^^ + ( n - c ■ (6-1) 

The field equations for ip a can be found by varying (6.1) with respect to ip a . We can then 
separate variables with the aid of similar expressions to those used in chapter 3 



J^a(j,^ e il0+i{wt+kz) _ j^ot* ^y^ e —il6—i(wt+kz) 



\ 

(6.2) 



/ 



q a are the charges of ip a with respect to the string generator. I can take any value which 
gives a single valued ip a , so it must be an integer if 2q a are odd and a half integer if 2q a 
are even. Using the gauge invariance of (6.1), we deduce that the angular dependence of 
the mass term is 

M ap (r,9) =C aP {r)e i{qa+q P )e (no summation) . (6.3) 
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We will define the following operators 



V y T = d r + — + diag h eA e , h eA e , eA e , eA e , 

2r V r r r r 
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(6.4) 



Using (6.2) the field equations reduce to 



(-V + V T + M) A = , (6.5) 

where (PyA)" = T>^A a . These have equal numbers of real and imaginary solutions. If 
/ = (which is only possible if 2q a are odd) these equations are reducible. We will consider 
real and imaginary solutions separately, and set A^ = ±Af and A" = ±^3- Defining 
A a = {Af,Af) T , (6.5) becomes 



with 



(-V + V T ± M) A = , 



V p= dr + l + l diag (-^ + eAe ,^-eA e ) . 
2r r \ r r 



(6.6) 



V = 



w + k 



( 



, M af3 = 



V 
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-1 



a 



(a/3) 



(6.7) 



w — k 

The upper and lower choices of sign correspond to real and imaginary solutions respectively. 
Now consider A T V 2 A. Using (6.5) and the definition of V, we find 



A T V 2 A = (w 2 - k 2 )\A\ 2 = A T {V T + M) 2 A . 
We will now define the constant orthogonal matrix 



(6.8) 



K 



—I 

1 



(6.9) 



Any physical solutions must be normalisable. We will scale A so that J \ A\ 2 rdr = 1. Now 
(6.8) implies 



w 



-k 2 = J A T {V T + M)K T K{V T + M)Ardr = J \K{V T + M)A\ 2 rdr . (6.10) 
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We have integrated by parts, and used the fact that M and T>t are symmetric. We see 
that w 2 > k 2 , with equality occuring only when (T>t + M)A = 0. Similar arguments 
apply to V. Thus there are no space-like fermion currents, and the only light-like ones are 
generalisations of zero modes. This is in contrast to the case of bosonic currents, where 
space-like currents exist and are the most favoured [58]. 

6.3 Bound States in the Abelian String Model 

For simplicity we will consider currents in the U(l) model described in chapter 1. The 
fermionic part of its Lagrangian is given by (1.16), which leads to the field equations (1.17). 
In this case there is only one fermion field, so the index a in (6.2) can be dropped. For a 
string with winding number n, q = n/2 and M = m{f(r)e m9 . 

Since there are no space-like currents we can set k = without loss of generality. By 
changing the signs of A3 and A4, w > solutions can be changed into w < solutions, thus 
only positive w need be considered. We will take the string's winding number (n) to be 1, 
so q = 1/2 in (6.2). By interchanging Ai^ and Ai^ negative / solutions can be obtained 
from the positive / solutions. Thus it is only necessary to look for / > solutions. 

We will use a variation of the shooting method to determine the values of w which 
have normalisable fermion solutions on the string. At large r the solutions of (6.5) have 
exponential behaviour. Two of them decay and so are acceptable. In the case of the small 
r solutions, only two of them give a normalisable state. 

Each of these 4 solutions can be numerically extended to some intermediate value of r 
(of order the string width). We can then see if any non-trivial combinations of the large 
and small r solutions match up there. Since the equations are linear, this is sufficient to 
determine if there are normalisable solutions for a given value of w. Thus we need only 
consider variations of w. The special / = case can be treated similarly. There is then just 
one large r and one small r acceptable solution. 

Figures 6.1 and 6.2 show the variation of the number of bound states with respect to the 
Yukawa coupling (mf/m s ). Each line corresponds to one real and one imaginary solution. 
For simplicity we have taken the Higgs and gauge field masses (m s and m v ) to be equal. 
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Figure 6.1: Spectrum of I = ... 3 fermion bound states. 

The string field profiles are shown in figure 1.2. 

Plots of the solutions for mf/m s = 2.3 are shown in figure 6.3. As expected they decay 
outside the string (which is approximately of radius 3 on the plots). 

6.4 Bound States in Other Theories 

We know that the right-handed neutrino zero modes on an 50(10) abelian string do not 
survive the electroweak phase transition (see section 3.3.1). Considering the form of the 
mass terms in the theory will give some insight into what happens to them. 

After the electroweak phase transition the left and right neutrino fields [y and v c ) receive 
masses from both <I>io and 3>i26- In the absence of cosmic strings, $io contributes a Dirac 
mass term m u = g E i]u and $126 contributes a Majorana mass term m G = g G rj G . Assuming 
m G 3> m u the neutrino mass eigenvalues are 

J nil + 4m. 2 , + m G 
m R = « m G 
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Figure 6.2: Spectrum of I = 4 ... 7 fermion bound states. 



rriL 



2 m G 



(6.11) 



The mass eigenstates are then approximately v c + eu and v — eis c , with e = m u /m G . 

Since |$i26| an d l^iol var y inside a cosmic string, the neutrino masses will too. In 
section 2.6.1 we showed that there is a region around the string core in which |3>io| is 
reduced but | ^>i26 1 takes its usual VEV. Thus mi will be lower in this region. At the centre 
of the string <j? 12 6 = and l^iol ~ m u / m G ( see figure 2.4), thus = ~ rr? u jra G there. 

We see that there is a potential well inside the string, which suggests that neutrino 
bound states will be present. As with the abelian string, we can investigate the existence of 
bound states by examining the approximate large and small r solutions of (6.5), and then 
trying to match them at intermediate r. The 8 small r solutions have the same behaviour 
as the small r zero mode solutions given by (3.12,3.13) and their right moving equivalents. 
Four of them have acceptable small r behaviour. Of the 8 large r solutions 4 decay outside 
the string (as required) if w 2 — k 2 < m|, so bound states could exist. If w 2 — k 2 > m 2 L no 
more than 2 of the solutions are acceptable, and so there are no bound states. Thus any 
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fermion solutions which are localised to the string will have w 2 — k 2 < m 2 L . The maximal 
current in this case will be tiny. 

At the phase transition the right-handed neutrino zero modes will mix with the left- 
handed neutrinos to produce the bound states. Any currents with rest mass greater than the 
left-handed neutrino mass will be able to escape from the string. The remaining current 
will be far to small to stabilise vortons, so they will certainly collapse. Furthermore we 
expect the bound state to be spread over the region of electroweak symmetry restoration. 
This is far larger than the size of a vorton, which is a couple of orders of magnitude greater 
than the GUT string radius [24]. The current on one part of the string will then interact 
with current on the opposite side of the loop, increasing the vorton's instability. Thus 
we conjecture that vortons in 50(10) will decay at the electroweak phase transition, even 
though currents may continue to exist. 

The electron and quark masses are also reduced inside the string, so they may have 
bound state solutions too. Their off-string mass is far greater than the neutrinos, so these 
currents may become large enough to have detectable effects. Since they are charged they 
will produce a wider range of effects than the neutrino currents. Even so, they will still be far 
smaller than the GUT currents that were present before the electroweak phase transition. 

The situation is likely to be different in super symmetric theories. Here, SUSY-breaking 
destroys the zero modes by mixing them with other zero modes, rather than massless 
fields. Again we expect the current to be limited by the off-string mass of the particles it is 
composed of. SUSY-breaking does not significantly alter this, so the current will not fly off 
the string as it could in the 50(10) model. Additionally, extra symmetry restoration does 
not occur, so the width of the string will not increase and the current on a vorton will not 
scatter off itself. This suggests vortons arising in supersymmetric theories are more stable 
than those in ordinary GUTs. Since the mass of the SUSY bound states will be of order 
the SUSY breaking scale, it is still possible the currents could decay into Standard Model 
particles. This will depend on the couplings of the theory. Hence fermion currents may 
provide constraints on SUSY models after all. 
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6.5 Summary 

The existence and form of fermion bound states and their corresponding currents on cosmic 
strings were investigated in this chapter. We have shown that only time-like and light-like 
currents can occur. Using numerical methods, the discrete spectrum of states for the abelian 
string model was determined. We found that it varied with the magnitude of the Yukawa 
coupling of the theory. For very low values there is just one bound state and a zero mode. 

Since the bound states can carry angular momentum, they will contribute to vorton 
formation. However, unless interactions between the states and the electroweak sector are 
suppressed, they will decay. 

We also speculated about the existence of massive fermion currents in models whose 
fermion zero modes are destroyed by a phase transition or supersymmetry breaking. We 
conjecture that when the zero mode mixes with a massless field, the size of the current 
will be reduced. If the energy scale at which this mixing occurs is much lower than that at 
which the string formed (as in the GUT considered), the current will be drastically reduced. 
Any vortons will then certainly decay. If the zero mode mixes with another zero mode (as 
in the SUSY theories considered), the current will not be significantly altered. Any such 
vortons must decay by another mechanism. 
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Chapter 7 



Summary and Future Work 

7.1 Summary 

In this thesis we have investigated cosmic strings in realistic particle physics theories. We 
have seen that these strings have a far richer microstructure than their counterparts in the 
simple abelian string model. This microstructure can lead to observational consequences, 
and will change the physical predictions of the theory. The additional properties of these 
strings can help to provide explanations for cosmological phenomena. They may also be 
used to constrain the underlying theory. 

We began by looking at cosmic strings in theories with several phase transitions. In 
chapter 2 we showed that a string formed at one phase transition could affect future phase 
transitions, and force their Higgs fields to take stringdike solutions. Even if stringdike 
solutions formed at such transitions would normally decay, the presence of the GUT de- 
fect stabilises them. The existence of these stringdike solutions leads to extra symmetry 
restoration. 

This effect was considered in detail for a realistic grand unified theory with an SO (10) 
symmetry group. After the electroweak phase transition there are a total of five distinct 
types of cosmic string in this model. Three of them restored electroweak symmetry in a 
region around the string core. The size of this region is of the same order as the electroweak 
scale, so it is far greater (about 10 14 times bigger) than the string width. Profiles of the 
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electroweak fields were found numerically. 

Previous work suggests that if fermions couple to the Higgs field of an abelian cosmic 
string, fermion zero modes (and hence massless currents) will exist on the strings [26]. 
In chapter 3 the existence of these currents in a general theory was investigated. We 
considered theories with more exotic cosmic strings whose Higgs fields can have more than 
one winding number. An index theorem giving the number and type of massless fermion 
currents was derived. The theorem depends only on the type of Yukawa terms present, 
and their angular dependence. It does not depend on the size of the Yukawa couplings, 
or the radial dependence of the string solution. This allows a wide range of theories to be 
examined with relative ease. We showed that massless currents can be created or destroyed 
at phase transitions. This provides a solution to the potential vorton problem. 

The index theorem was applied to the strings discussed in chapter 2. We found that 
for topologically stable strings, only the abelian ones had massless currents at high tem- 
peratures. In this case the current carriers were right handed neutrinos. Following the 
electroweak phase transition, no topologically stable strings admitted massless currents. 
Thus fermion superconductivity is unlikely to occur in this model at low temperatures. 
The model provides an example of massless current destruction. Some of the strings with 
higher winding numbers did have fermion currents, although such strings are unlikely to be 
stable. 

The idea that particle physics theories are supersymmetric above the electroweak scale 
has become increasingly favoured in recent years. In chapter 4 the implications of super- 
symmetry for abelian cosmic strings were investigated. Two methods of gauge symmetry 
breaking were considered, giving two types of string. The strings broke (fully in one theory, 
partially in the other) supersymmetry in their cores. Applying supersymmetry transforma- 
tions did not therefore leave the solution invariant. However, since the field equations were 
still invariant under supersymmetry, the transformed solutions still solved them. These new 
solutions corresponded to strings with fermion zero modes on them. Furthermore these so- 
lutions were expressed in terms of the background string fields, whose properties have been 
extensively studied. 
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Supersymmetry is clearly broken in the present universe. We examined the effects of soft 
supersymmetry breaking terms on cosmic strings in chapter 5. A theory with abelian and 
nonabelian cosmic strings was considered. Soft terms are those which break supersymmetry, 
but do not give quadratically divergent contributions to the electroweak Higgs field mass. 
Some of them destroyed the zero modes. There is some similarity between this mechanism 
and the destruction of high temperature neutrino zero modes in chapter 3. The effects of 
soft supersymmetry breaking on the abelian strings in chapter 4 were also considered. The 
results were the same for the F-term symmetry breaking, but the zero modes on strings 
formed by D-term symmetry breaking were unaffected. 

Lastly, in chapter 6, we investigated the spectrum of massive fermion bound states and 
currents in an abelian string model. We found that unlike the massless currents, the number 
of massive currents is dependent on the strength of the couplings in the theory. We also 
showed that space-like fermion currents do not exist on cosmic strings in any model. 

We have shown that cosmic strings in realistic models can be significantly different to 
those in simple abelian models. If quantum field theory is supersymmetric at high temper- 
atures cosmic strings are certain to have conserved currents. This can radically alter the 
cosmological implications of the strings. We have also found that the properties of realistic 
strings are not constant, but change as the universe passes through phase transitions (or 
as supersymmetry is broken). For example, the stability of currents on the strings can be 
changed. 

7.2 Future Work 

There are clearly many extensions of the work on fermion bound states started in chapter 6. 
The destruction of massless fermion currents at gauge and supersymmetry breakings has 
been seen in chapters 2 and 5. It seems likely that they turn into bound states, which may 
then decay. I will examine this conjecture in more detail. If such currents do decay, analysis 
of the bound states will help determine the decay products. Such decays may help provide 
a mechanism for baryogenesis. 

Although the existence of conserved fermion currents on strings could stabilise loops, 
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there is still debate about this. The currents are not topologically conserved, and there are 
possible mechanisms that could destabilise them. One possibility is the interaction of the 
current carriers with plasma particles or other current carriers. These interactions can also 
create currents. The presence of fermion bound states on a string may alter these processes. 
I will investigate these ideas. 

The only GUTs considered in chapter 2 involved SO (10) (or a subgroup), but many of 
the results can be generalised to other theories, such as Chapter 3 could be extended 
in a similar way. I am currently investigating the zero modes in an axion cosmic string 
theory [59], which can arise naturally in superstring theory. It is possible that phase 
transitions will also affect gauge boson currents, in a similar way to fermion currents. Such 
currents are a common feature of nonabelian theories. 

Another possible area of investigation is the effect of the other symmetry restora- 
tions (particularly the SU(5) restoration) on monopoles. In the S*O(10) GUT considered, 
monopoles form when SU(5) is broken [60]. It is credible that it would be energetically 
favourable for monopoles to sit on nonabelian strings, since this would reduce the variation 
in $45, and so possibly the total energy. If this does happen, the likelihood of monopole 
collision will be greatly increased, as collisions in one dimension are far more frequent than 
in three. Since they are Z2 monopoles (and hence monopoles are topologically equivalent 
to anti- monopoles), this results in a higher annihilation rate. While monopoles occur in 
most GUTs, their properties conflict with observations, so a mechanism is required to get 
rid of them. Since cosmic strings could reduce the number of monopoles they could solve, 
or help to solve, this problem. I will see if this mechanism works, and then determine its 
effectiveness. I will also see how the presence of a gauge boson current on the cosmic string 
affects it. This mechanism has similarities to the Langacker-Pi mechanism [61], in which 
strings form that link the monopoles together. The monopoles are then ends to the strings. 
These rapidly contract, pulling the monopoles into each other, causing them to annihilate. 
This mechanism is different because the strings are attached to the monopoles for topolog- 
ical reasons rather than just dynamical ones. While the Langacker-Pi mechanism is very 
efficient, it requires the breaking and later restoration of electromagnetism, which is hard 
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to reconcile with experiment. 

One more aspect of GUT strings that I intend to examine is the interaction of several 
strings. Collisions between abelian strings have been extensively studied, while those be- 
tween the more exotic string solutions have received less attention. I will look at these in 
more detail, particularly in realistic GUTs. Previously, strings that form in a G — ► H x D 
symmetry breaking, where the discrete group D is nonabelian, have been considered. These 
strings cannot intercommute for topological reasons [62]. This leads to a very different evo- 
lution of cosmic string networks [63]. Unfortunately such strings are generally Alice strings, 
and so are not physical. It is possible similar effects will occur in more realistic theories, 
although for dynamical rather than topological reasons. Even the simplest abelian strings 
do not intercommute in some cases [29]. If strings do not intercommute, loop production 
will be suppressed. A different mechanism for energy loss is then needed if the universe 
is not to become string dominated. If there are significant differences with GUT strings, 
they will have implications for the evolution of string networks. This is particularly im- 
portant given that recent simulations with abelian strings [64] indicate that they produce 
too little temperature anisotropy at small angular scales to explain the observed data. The 
cosmology of more exotic cosmic strings remains an open question. 

The effects of soft SUSY breaking have only been partially explored in chapter 5, and I 
will be considering them further. It is not clear what form the minimum of a general scalar 
potential takes when soft SUSY breaking terms are present. It appears that it may have 
a flat direction (in addition to the one arising from the gauge symmetry). This could be 
of interest, since it has been suggested that it could lead to a force between strings. Flat 
directions are not uncommon in SUSY theories, and a more approachable theory which also 
has strings has been considered in [65]. If a similar mechanism can be found for monopoles, 
it could provide a possible solution to the monopole problem. It is also possible the flat 
direction will give some kind of inflation (rapid expansion of the universe, which could solve 
many cosmological problems), as it does in unbroken supersymmetry. I hope to extend this 
work to supergravity theories. 

Another interesting property of SUSY is that inflation and cosmic strings frequently 
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occur together. When cosmic strings are formed because of F-terms in the potential, 
hybrid inflation also occurs. When a D-term is used to produce inflation, the breaking of 
the U(l) will result in string formation. I will look at the implications these two theories 
have for each other. In particular, the cosmological predictions are likely to be much richer 
than those resulting from strings or inflation alone. 
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Appendix A 



An SO (10) Grand Unified Theory 

A.l Fermion and Gauge Fields 

Under 50(10), all left-handed fermions transform under one representation, and right- 
handed fermions transform under its conjugate [13]. It is convienient to use just one repre- 
sentation. This can be achieved by using the charge conjugates of the fermions (tp R = O^, 
i/j c l = C$1). The charge conjugate of a right handed fermion transforms as a left handed 
fermion, and vice versa. Thus *f>L = ipL + $% is left handed. For 50(10) this definition 
is necessary, as well as convienient. could be gauge transformed to ip R , so any gauge 
invariant quantities will have to involve just ^ l and * r (right-handed equivalent of ^l)- 
However, is superfluous, since it is equal to ia 2 ^* L = ^> C L , so the theory can be described 
entirely in terms of For the electron family, it can be written as 

= ( u i^ u 2,u 3 ,u e ,d 1 ,d 2 ,d 3 ,e^ ,dl,d9,,dl,e + ,-u c 1 ,-U2,-ul,-u^) T , (A.l) 

where di = dn, d\ = ia 2 d* R , etc. so all the fields are left handed. The other two families of 
fermions can be described similarly. Adapting work by Rajpoot [66], the gauge fields can 
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be expressed explicitly as 16 x 16 matrices, which act on the fermion fields. 
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Xci Xr 



X S2 
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53 



53 
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(A.3) 



(A.4) 



L 51 ^ v 52 

G is a 3 x 3 matrix of containing the gluon fields. It is hermitian, and so H is too. The 
other fields are contained in the diagonal matrix Ai6 



Ai6 = diag((^ + Wl) o , -3& + W L \ " Wl)^ -3 



\V6 



'76 



wl, 



(-S + ^)3' 3 * + ^(-76-^) 3 ' 3 S 

diag((s + 2a + 2z)s, —3s + 4z, (s — 3z — 0)3, —3s — 3a - 
(—3s + a — z) 3 , s + 3a — 3z, (s — 2a + 2z) 3 , 5s) 



wl 



z, 



(A.5) 



The subscripts indicate repeated values, and 

1 



-Wl + f-B> 



(A.6) 
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B = Vf^i + Vf^', (A.7) 
z = \(wt- f-B^ , (A.8) 

a = I ( Wl * + \fl B ) ■ (A - 9) 

Z = vTOz and A = \/6a are the unrenormalised electroweak Z° boson and photon respec- 
tively. Z' = —y/Ws is a high energy 50(10) boson. The generator P is obtained by putting 
s = 1 and a = z = in the expression for Ai6. The substitutions s = z = and a = 1/3 
give the charge operator. 



A. 2 Higgs Fields 

The electroweak Higgs field transforms under the 5_2 and 52 representations of SU(5). 
They are contained in the 10 of 50(10), and so the components of <3?io can be expressed as 
symmetric products of spinors transforming under the 16 representation. Thus 3>io can be 
expressed as (ffiHj- + (f>Z^u > where H% (a = 0, ±, 1, 2, 3), are the five components of 5_2, 
and H" are the corresponding components of 82- (H®,H£) and (H®,H~) form SU(2)l 
doublets, while H l d and H % u (i = 1,2,3) form an SU(3) C triplet and anti-triplet. 
Expressing these components of 10 in terms of symmetric products of 16s gives 

Hi = iKujXu^s + iOx^s}, 

H~ = -kiiVj * dj) s + (e~ xu c ) s ] , 

K = -\hjk(uj x d k ) s + (d c i x v c ) s - (u c i x e + ) s ] , 

H° d = iKdjxtys + ie+xns] , 
H+ = i[(^xdp 5 + (e + x%] , 

H d = lhjk(u c j x d c k ) s + (di x % - (Hi x e~) s ] , (A.10) 
where v is the basis vector corresponding the v field, etc. The gauge fields can also be 
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expressed as 10 x 10 matrices, which will that act on (^,0"). 
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~ X S2 


X S1 





X '3 + 


-yf 


-X? 


~ X 2 + 


- x 's + 
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-2^) , Wl + Wl -Wt + Wl 



(A.11) 



(A.12) 



(A.13) 



-w£-w&,w£-m 



= diag ((—2s — a + z)%, —2s + 3a + z, —2s — 4z, 

(2s + a - z) 3 , 2s - 3a - z, 2s + 4z) 



(A.14) 



The second Higgs field $45 is in the 24q component of the 45 representation, and its usual 
vacuum expectation value is proportional to the generator of the B field. In the absence of 
a string, <J> 12 6 is proportional to (ei x ei)s, with ei = — v c . 



A. 3 Fermion Masses 

The masses of the fermions arise from Yukawa couplings to the Higgs fields. These must of 
course be Lorentz and gauge invariant. Since the theory contains only left handed spinors 
the only possible Lorentz invariant mass terms are Majorana masses (tplipL + V'lV'l + 
(right-hand terms)). The Majorana masses transform as a product of 16s, and so can be 

113 



coupled to similarly transforming Higgs fields, <E>i26 and 5>io, but not $45. This gives the 
fermionic Lagrangian (3.36). We have considered only one family of fermions for simplicity. 

Of course, since (3.36) is invariant under 50(10), it must be invariant under SU(5) 
as well. Thus ^126 can on b/ couple to SU(5) singlets (i.e. products of the conjugate 
neutrino field). The two components of QJff couple to 5_2 and 82 products of fermions. 
Under 577(5), the remaining fermions transform under 10i and 5_3 representations. To 
find allowable mass terms, products of these representations need to be expressed in terms 
of irreducible representations. 10i x 10i and I5 x 5_3 both contain S2S, and 10i x 5_3 
contains a 5_2- Thus, for the usual VEVs of the Higgs fields, the mass terms are written 
in terms of particle fields as 

V L H° d if c L = § [d\ia 2 dT + e-Ua 2 e + * + df ia 2 d* + e + Ua 2 e~* 

= \ diLdiR + e~ L e~ R , (A.15) 



and similarly for H®, so 



^L^lon = it d tL d tR + ele R + ^[u iL u iR + v L v R \ , (A.16) 



4 



4 



^l^I26^ C l = r) G v c Ha 2 v c * = r )G v T R io 2 v R . (A.17) 

This model, unlike the standard model, has non-zero neutrino masses. If the v^i^vr term 
is much larger than the v^vr term, the mass eigenstates will be approximately vl and vr, 
and have very small and very large mass eigenvalues respectively, giving an almost massless 
left-handed neutrino. 
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